Real-time noise-suppression method and apparatus therefor

ABSTRACT

A thorough error suppression signal measurement system ( 20 ) having a transmitter ( 300 ) for propagating a transmission signal to a first probe  100 , through a device under test ( 26 ), and into a second probe ( 200 ), and for propagating reference signals to the probes ( 100,200 ). The probes ( 100,200 ) extract normalization signals from the reference signals, exchange specific ones of the normalization signals, and combine the normalization signals with data signals derived from the transmission signal to form receiver signals. The probes ( 100,200 ) propagate the receiver signals to a receiver ( 400 ), where the signals are gain-ranged, digitized, normalized, and compensated for phase-noise.

RELATED INVENTIONS

The present invention claims priority under 35 U.S.C. §119(e) to: “An Accurate, Real-Time, Total-Error-Suppression Technique (ARTTEST Method) and Its Application to High-Accuracy, Network-Analyzer Calibration and Measurements,” Provisional U.S. Patent Application Ser. No. 60/101,225, filed 21 Sep. 1998, which is incorporated by reference herein.

The present invention is a divisional (DIV) of “Real-Time Error-Suppression Method And Apparatus Therefor,” U.S. patent application Ser. No. 09/400,220, filed 21 Sep. 1999, which is incorporated by reference herein.

TECHNICAL FIELD OF THE INVENTION

The present invention relates to the field of signal measurement. More specifically, the present invention relates to the field of integral and simultaneous signal measurement and measurement device calibration.

BACKGROUND OF THE INVENTION

For accurate instrumentation, it is desirable to fully understand the characteristics of the devices used. With any given one of these devices, i.e., a “device-under-test” or DUT, the specific characteristics of that DUT are not fully understood. As a simplistic example, a given “47 kΩ” resistor would rarely have a value of exactly 47,000Ω. To properly understand the operation of a circuit containing such a resistor, a knowledge of the actual resistance (e.g., 46,985.42Ω) would be helpful. It should be understood that a DUT may encompass a wide range of electrical components, equipment, and systems. A typical DUT may be a filter, an amplifier, a transmitter, a receiver, or any component, group of components, circuit, module, device, system, etc.

The drift of components, filters, amplifiers, and other signal-conditioning circuitry typically limits the accuracy of electronic measurements. A measurement system may advertise a large dynamic range and very-high resolution. However, the full dynamic range and resolution may not be realizable because of errors inherent in the system. Current measurement systems have not been demonstrated to have verifiable, full-scale uncertainties of better than 0.1 percent over all types of errors.

Real-world measurement instruments tend not to be perfect. This imperfection will affect the accuracy of the resulting measurements. This accuracy is dependent upon measurement errors. Such errors may be classed as static (systematic) errors and dynamic errors.

Static errors are repeatable, time-invariant system errors. That is, static errors do not vary over time. Static errors result from the nonideal aspects of a system. These errors are repeatable as long as no changes are made to the system. Static errors include directivity errors, source-mismatch errors, load-mismatch errors, reflection and transmission tracking errors, isolation or cross-talk errors, etc.

Static errors may be reduced through the use of precision components and circuits. However, no matter how precisely a circuit is designed, there will still be some level of static error present. Since static errors are repeatable, they can be suppressed using various static error suppression techniques, such as twelve-term error modeling, known to those skilled in the art. Twelve-term error modeling, typically employed with standard network analyzers, can account for directivity errors, source-mismatch errors, load-mismatch errors, tracking errors, and isolation errors.

Before error modeling may be employed, the error coefficients of the requisite equations must be calculated by making a set of measurements on a set of known loads meeting precise standards. A sufficient number of precise standards must be used in order to determine the various error coefficients in the error model. A common static-error calibration technique is the Short, Open, Load, and Through (SOLT) technique. The SOLT calibration technique yields better than 0.1 percent accuracy for static errors. However, dynamic errors limit the actual accuracy to less than this. An alternative static-error calibration technique, also well known to those skilled in the art, is the Through, Reflection, and Load (TRL) technique. The TRL calibration technique yields significantly better static-error accuracy than the SOLT technique at the cost of calibration complexity. There are a number of other well-known static-error calibration techniques that may be used to advantage in specific instances.

However, removing a significant amount of static error achieves little if other types of errors are not also addressed. Dynamic errors are time-varying errors. That is, dynamic errors change over time, often in an unpredictable manner. Such errors may be attributable to a number of different sources. For example, component-drift, physical-device (e.g., cables, connectors, etc.) errors, phase-noise errors, random noise, etc.

Measurement systems exhibit several types of dynamic errors. Phase-noise errors and random errors, while inherently dynamic (i.e., time-variant) are special cases independently discussed hereinafter.

Component-drift errors may be either long-term or short-term. Long-term component-drift errors, however, are typically due to aging of the components, with resultant variations in component specifications.

One type of short-term component-drift error, source drift error, is usually attributable to thermal or mechanical variations within the system, and may include both amplitude and phase fluctuations of the output wave.

Another type of short-term component-drift error, receiver drift error, is associated with a data receiver. This may be due to drift in components such as amplifiers, filters, and analog-to-digital converters. Receiver-drift error may also appear as time-varying gain and phase drift in the received signals.

Thermal variations may also lead to physical expansion of passive components within the system. At high frequencies, such expansion may lead to appreciable phase errors. In applications where the DUT is located at a considerable distance from a transmitter and/or receiver, there may be a number of time-varying errors associated with the connections between the DUT and the transmitter and receiver. These may include amplitude and phase-drift errors in the amplifiers or errors associated with the modulation and demodulation circuitry. Systems in which such errors become significant include systems utilizing propagation media other than traditional cables (e.g., the atmosphere, space, the earth, railroad tracks, power transmission lines, the oceans, etc.).

Dynamic physical errors result from physical changes in the test setup. One example of a physical error is connector repeatability. As one connects and disconnects the DUT, there will be reflection and transmission errors associated with any nonrepeatability of the connectors. The severity of the connector repeatability error is related to the type of connector, the condition of the connector, and the care with which the user makes the connection.

Cable-flexure errors are another type of dynamic physical error. Cable-flexure errors appear as one moves the cables to connect or disconnect a DUT or perform a calibration. Time-variant phase errors associated with the relaxation of the cable can occur for a period of time after the cable has been flexed.

Phase noise (jitter) is directly related to the frequency stability of a signal source. In a perfect sinusoidal oscillator, all the energy would lie at a single frequency. Since oscillators are not perfect, however, the energy will be spread slightly in frequency. This results in a pedestal effect. This effect, referred to as phase noise, is more severe at higher frequencies. Phase noise is a performance-limiting factor in applications where a weak signal must be detected in the presence of a stronger, interfering signal.

Random or white noise is common in measurement systems. Random noise includes thermal noise, shot noise, and electromagnetic interference. Random noise may appear as random data errors. Traditional and well-known techniques of random error suppression utilize some form of oversampling to determine the correct data and suppress the random errors.

Calibration frequency is also a problem in conventional signal measurement systems. Typically, high-accuracy measurement systems employing manual calibration are calibrated periodically. The interval between calibrations may be hourly, daily, weekly, monthly, quarterly, or even yearly. This technique produces a steadily decreasing accuracy that progresses over the inter-calibration interval. Additionally, drift errors occurring during the inter-calibration interval are uncompensated. Such drift errors tend to accumulate. Therefore, measurements taken shortly before calibration may be suspect. Exactly how suspect such measurements may be depends upon the length of the inter-calibration interval and the amount of drift involved.

Many state-of-the-art measurement systems employ automatic-calibration techniques. Some automatic-calibration systems calibrate at specified intervals. Such systems suffer the same decreasing accuracy as manual-calibration systems.

Other automatic-calibration systems calibrate at the beginning and end of each measurement cycle. The use of frequent nonsimultaneous calibration procedures does increase overall accuracy, but may also greatly increase the cost of measurements and prevents measurements while the frequent calibration procedures are taking place.

All calibration procedures discussed above are nonsimultaneous. That is, the calibration procedures do not occur simultaneously with measurement. Nonsimultaneous calibration procedures are incapable of correcting or compensating for dynamic errors, such as component drift, occurring during measurement.

Current technology demands increasingly small operational errors. High accuracy is therefore a growing need of instrumentation users. Measurement systems utilizing simultaneous calibration are useful for applications requiring high-accuracy measurements. That is, systems are needed that calibrate themselves and measure data simultaneously. Such systems are said to employ dynamic error suppression. That is, such systems are able to compensate for dynamic (time-variant) errors by continuously calibrating themselves while simultaneously performing the requisite measurements.

As current technology drives operational frequencies higher and higher, phase noise (i.e., signal jitter) increases in importance. A definite need exists, therefore, for systems employing phase-noise error suppression. That is, for systems employing some means of compensating for signal jitter. This is especially important in polyphase-constellation communications systems where phase noise may cause misinterpretation of the signal phase point (i.e., the data) at any given instant.

Measurement systems for state-of-the-art technology also desirably compensate for random errors.

SUMMARY OF THE INVENTION

Accordingly, it is an advantage of the present invention that a real-time error-suppression technique is provided.

It is another advantage of the present invention that an apparatus for real-time error-suppression is provided.

It is another advantage of the present invention that a technique is provided to effect thorough error suppression, i.e., static error suppression, dynamic error suppression, phase-noise error suppression, and random error suppression.

It is another advantage of the present invention that simultaneous calibration is provided to allow dynamic error suppression.

It is another advantage of the present invention that signal normalization is provided to effect thorough error suppression.

It is another advantage of the present invention that simultaneous measurement of both data and normalization signals is provided to allow thorough error suppression.

It is another advantage of the present invention that a phase-noise reference signal is provided to effect phase noise suppression.

It is another advantage of the present invention that frequency offsetting is provided to effect discrete signal normalization.

It is another advantage of the present invention that the normalization signals are created from the data signal.

It is another advantage of the present invention that a bidirectional link is provided between input and output transmission signals to effect system calibration.

The above and other advantages of the present invention are carried out in one form by a signal measurement method incorporating real-time error suppression, which method includes combining a first signal and a second signal offset in frequency from the first signal to produce a third signal, propagating the third signal over a signal path capable of inducing errors to produce a fourth signal, separating the fourth signal into a fifth signal and a sixth signal, comparing the sixth signal against the second signal to determine a correction factor, and applying the correction factor to the fifth signal to extract a replica of the first signal therefrom.

The above and other advantages of the present invention are carried out in another form by a signal measurement system incorporating real-time error suppression, which system incorporates a signal source generating a first signal, a first mixer offsetting the first signal to produce a second signal, a second mixer offsetting the first signal to produce a third signal, a combiner combining the second and third signals to produce a fourth signal, a propagation medium propagating the fourth signal over a potentially corruptive signal path to produce an analog fifth signal, an analog-to-digital converter converting the analog fifth signal into a digital sixth signal, and a processor containing a separation routine separating the sixth signal into seventh and eighth signals, a comparison routine comparing the eighth signal against the third signal to determine a correction factor, and a correction routine applying the correction factor to the seventh signal to extract a replica of the second signal therefrom.

The above and other advantages of the present invention are carried out in another form by a signal measurement method for a device-under-test having a first port and a second port, which method incorporates coupling a first probe to the first port, coupling a second probe to the second port, producing a first signal, extracting a second signal from the first signal, propagating the second signal into the first port, deriving a third signal from the second signal within the device-under-test, propagating the third signal from the second port, deriving a fourth signal from the second signal, deriving a fifth signal from the first signal, deriving a sixth signal from the third signal, deriving a seventh signal from the first signal, reciprocally propagating the fifth and seventh signals from the first and second probes to the second and first probes, respectively, over a common signal path, combining the fourth and seventh signals to produce an eighth signal, and combining the sixth and fifth signals to produce a ninth signal.

The above and other advantages of the present invention are carried out in another form by a signal measurement system incorporating a device-under-test having a first port and a second port, a first probe coupled to the first port, a second probe coupled to the second port, a signal source generating a first signal, an input circuit extracting a second signal from the first signal and propagating the second signal into the first port, an output circuit deriving a third signal from the second signal and propagating the third signal from the second port, a first mixer producing a fourth signal from the second signal, a second mixer producing a fifth signal from the first signal, a third mixer producing a sixth signal from the third signal, a fourth mixer producing a seventh signal from the first signal, a trans-probe propagation medium reciprocally propagating the fifth and seventh signals from the first and second probes to the second and first probes, respectively, over a common signal path, a first combiner combining the fourth and seventh signals into an eighth signal, and a second combiner combining the sixth and fifth signals into a ninth signal.

BRIEF DESCRIPTION OF THE DRAWINGS

A more complete understanding of the present invention may be derived by referring to the detailed description and claims when considered in connection with the Figures, wherein like reference numbers refer to similar items throughout the Figures, and:

FIG. 1 shows a block diagram of a signal measurement system realized as a network analyzer in accordance with a preferred embodiment of the present invention;

FIG. 2 shows a simplified schematic diagram of a transmitter of the network analyzer of FIG. 1 in accordance with a preferred embodiment of the present invention;

FIG. 3 shows a simplified schematic diagram of a first probe of the network analyzer of FIG. 1 in accordance with a preferred embodiment of the present invention;

FIG. 4 shows a simplified schematic diagram of a second probe of the network analyzer of FIG. 1 in accordance with a preferred embodiment of the present invention;

FIG. 5 shows a simplified schematic diagram of a receiver of the network analyzer of FIG. 1 in accordance with a preferred embodiment of the present invention;

FIG. 6 shows a simplified schematic diagram of a generic mixing circuit in accordance with a preferred embodiment of the present invention;

FIG. 7 shows a simplified schematic diagram of a processor of the network analyzer of FIG. 1 in accordance with a preferred embodiment of the present invention;

FIG. 8 shows a signal-flow graph of the signal-flow path of a signal of the network analyzer of FIG. 1 in accordance with a preferred embodiment of the present invention;

FIG. 9 shows a general signal-flow graph in accordance with a preferred embodiment of the present invention;

FIG. 10 shows a two-port 12-term error-correction model in accordance with a preferred embodiment of the present invention;

FIG. 11 shows a simplified schematic diagram of a reciprocal amplification circuit for a Link channel of the network analyzer of FIG. 1 in accordance with an alternative preferred embodiment of the present invention;

FIG. 12 shows a chart depicting relative frequencies for signals in the reciprocal amplification circuit of FIG. 11 in accordance with an alternative preferred embodiment of the present invention;

FIG. 13 shows a simplified schematic diagram of a circular mixer circuit for the network analyzer of FIG. 1 in accordance with an alternative preferred embodiment of the present invention;

FIG. 14 shows a simplified schematic diagram of a mix/demix circuit in accordance with an alternative embodiment of the present invention; and

FIG. 15 shows a simplified schematic diagram of an intermodulation distortion checking and removal circuit in accordance with a preferred embodiment of the present invention.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

An approach to signal measurement described herein is referred to below as a Thorough Error Suppression (TES) technique. TES applies a combination of techniques to a variety of significant and diverse errors typically encountered in signal-measurement systems. In one preferred embodiment, the TES technique combines static error suppression (SES), dynamic error suppression (DES), phase-noise error suppression (PNS), and random error suppression (RES).

The TES technique is generally applicable to a wide variety of measurement problems and can lead to measurements of greater accuracy than are possible using conventional techniques. The TES technique can be applied in many fields in which high accuracy measurement is limited by conventional measurement technology. These fields include communications, radar, detection, monitoring, imaging, etc.

FIG. 1 shows a block diagram of a signal measurement system 20 realized as a TES network analyzer 22 in accordance with one preferred embodiment of the present invention. TES network analyzer 22 has been chosen as an exemplary embodiment of signal measurement system 20 because TES network analyzer 22 incorporates a significant majority of the elements of the TES technique. Those skilled in the art will appreciate that many apparatuses other than a network analyzer may use some or all of the techniques discussed herein. The use of TES network analyzer 22 is an exemplary application, therefore, it neither indicates an application preference nor constitutes an application restriction in the use of either signal measurement system 20 or the TES technique.

The elements of signal measurement system 20 are shown in FIG. 1. System 20 incorporates a signal source 24, network analyzer 22, and a device-under-test (DUT) 26. Network analyzer 22 itself incorporates an analysis unit 28, a first probe 100, and a second probe 200. First and second probes 100 and 200 are configured to be coupled to a first port 30 and a second port 32 of DUT 26, respectively.

Analysis unit 28 contains a transmitter 300 configured to transmit various signals to first and second probes 100 and 200, hence to DUT 26, a receiver 400 configured to receive and process signals from first and second probes 100 and 200, and a processor 500 configured to analyze those received signals. Transmitter 300 is coupled to first and second probes 100 and 200 by high-frequency propagation media 102 and 202, and by low-frequency propagation media 104 and 204, respectively. High-frequency propagation media 102 and 202 merely convey higher frequency signals than do low-frequency propagation media 104 and 204, respectively. Those skilled in the art will appreciate that nothing requires high-frequency propagation media 102 and 202 to be different from low-frequency propagation media 104 and 204. The use of similar or substantially identical propagation media does not depart from the spirit of the present invention.

This discussion assigns reference numbers to components according to the following scheme. Reference numbers below 100 are assigned to components that are not a part of network analyzer 22. Reference numbers in the 100s are assigned to components that are a part of or directly associated with first probe 100. Reference numbers in the 200s are assigned to components that are a part of or directly associated with second probe 200. Reference numbers in the 300s are assigned to components that are a part of or directly associated with transmitter 300. Reference numbers in the 400s are assigned to components that are a part of or directly associated with receiver 400. Reference numbers in the 500s are assigned to components that are a part of or directly associated with processor 500. Reference numbers in the 600s are assigned to components that are a part of or directly associated with a reciprocal amplification circuit 600 in an optional embodiment of the present invention. Reference numbers in the 700s are assigned to components that are a part of or directly associated with a circular-mixer circuit 700 in an optional embodiment of the present invention. Reference numbers in the 800s are assigned to components that are a part of or directly associated with a mix/demix circuit 800 in an optional embodiment of the present invention. Reference numbers in the 900s are assigned to components that are a part of or directly associated with an intermodulation distortion checking and removal circuit 900 in an optional embodiment of the present invention.

In the preferred embodiment, the TES technique combines static, dynamic, phase-noise, and random error suppression. In the exemplary network-analyzer application, the design goal for TES network analyzer 22 was to achieve at least a 0.1 percent (−60 dB) measurement accuracy for DUT 26 with insertion and return losses varying over a dynamic range of 60 dB (i.e., a range of up to 120 dB). The main limiting factor in the achievable accuracy for TES network analyzer 22 resides in the absolute accuracy of the standards that are employed during static calibration. If a Short, Open, Load, and Through (SOLT) calibration technique is used, this accuracy would be approximately 0.1 percent. When a more accurate static calibration technique is being employed, such as the Through, Reflection, and Load (TRL) technique, or when only relative measurements are desired, the TES technique may be applied to achieve a measurement accuracy exceeding 0.1 percent.

TES network analyzer 22 uses signal-normalization techniques to achieve both static and dynamic error suppression for all components through which both data and normalization signals pass. TES network analyzer 22 is designed so that both data and normalization signals flow through all paths, e.g., transmitter 300, receiver 400, and link paths (to be discussed), with the exception of specific paths in first and second probes 100 and 200, the path through DUT 26, and errors caused by connector repeatability. This allows signal normalization techniques to be used to suppress static and dynamic errors in system 20.

The components in the above-mentioned specific paths in first and second probes 100 and 200 can be temperature stabilized so that they do not contribute to the system drift error. Furthermore, the static errors associated with the components in probes 100 and 200 are removed via 12-term calibration (discussed hereinafter). Since only the transmission signal sees DUT 26, the desired DUT characteristics are still present after normalization.

In order to minimize the nonremovable, time-varying errors associated with connector repeatability, it is desirable to use high-quality connectors, properly maintain the connectors, and to use a consistent, correct connection technique. When high-quality connectors are used appropriately in standard network analyzer measurements, connector-repeatability errors are reduced to below the level of residual systematic errors after correction.

The use of signal normalization for static and dynamic error suppression is illustrated by the following examples. Assume that two sinusoidal signals reside at different frequencies that are closely spaced relative to the component bandwidth. In a first case, with an assumption of no drift during the time of measurement, the two signals are simultaneously passed through a linear component whose properties vary between one measurement and the next. Both of these signals will then contain the same time-dependent component transfer function. If the composite time-domain signal is sampled by an A/D converter, processed using a numerical fast Fourier transform (FFT), and the resulting frequency-domain signals are extracted from the appropriate frequency bins, then the time-dependent transfer function can be removed via signal normalization, i.e., by computing the ratio of the two complex-valued, frequency-domain signals. Note that the normalization also removes the associated static errors.

In a second, more complex case, with an assumption of drift (i.e., amplitude or phase variations that occur in the component transfer function) during the time of measurement, the signals are again passed through the linear component. Drift appears as amplitude and phase modulation of the signals. As long as the rate of drift is not too large, the resultant signal spreading will not interfere with the neighboring signals. In such cases, both signals are affected the same by the amplitude and phase modulation, so static and dynamic error suppression can once again be accomplished by signal normalization. Signal normalization is executed in software in processor 500.

Since the incident, reflected, and transmitted signals are all to be measured in order to characterize a 2-port DUT in terms of its scattering-parameter matrices (S-matrices), TES network analyzer 22 includes a number of different normalization signals. Because TES network analyzer 22 has three primary components (first probe 100, second probe 200 and analysis unit 28) with a large number of signals distributed therein, this discussion utilizes a simple schema for signal nomenclature:

-   _(M)XX_(IJ)     where XX is the primary signal descriptor, e.g.,     -   ComF₁₁ port-1 fixed-common signal; and     -   Inc₂₃ port-2 composite incident signal;         where prefix subscript M is a secondary signal descriptor, e.g.,     -   _(R)Dat₁₁ port-1 reflected (raw) data signal;     -   _(O)Dat₁₁ port-1 offset data signal; and     -   _(M)Dat₁₃ port-1 measured data signal; and         where suffix subscript I is the signal reference and suffix         subscript J is the signal location, and where the values of I         and/or J may be:         -   1=first probe 100;         -   2=second probe 200;         -   3=analysis unit 28 or signal source 24;         -   P=first or second probe 100 or 200; and         -   J=P or analysis unit 28;             e.g.,     -   _(L)ComV₁₁ port-1 variable-common local oscillator signal in         first probe 100;     -   _(L)ComV₁₃ port-1 variable-common local-oscillator signal in         analysis unit 28;     -   Link_(1J) port-1 link signal in probe 100, probe 200, and/or         analysis unit 28; and     -   Link_(PJ) port-1 or port-2 link signal in probe 100, probe 200,         and/or analysis unit 28.

Signal source 24 generates or otherwise provides source signal _(S)Tx₃₃. Source signal _(S)Tx₃₃ is then processed by transmitter 300 within analysis unit 28 to produce either a port-1 transmission signal Tx₁₃ or a port-2 transmission signal Tx₂₃. Signal Tx₁₃ or Tx₂₃ is then transmitted from transmitter 300 to first or second probe 100 or 200 through high-frequency propagation media 102 or 202, respectively, where it arrives as a transmission signal Tx₁₁ or Tx₂₂.

Assuming forward signal measurement (i.e., the signal to be measured passes from first probe 100, into first port 30, through DUT 26, out of second port 32, and into second probe 200), then first probe 100 derives a forward DUT input signal _(I)Tx₁₁ from port-1 transmission signal Tx₁₁ and injects input signal _(I)Tx₁₁ into first port 30 of DUT 26. Within DUT 26, input signal _(I)Tx₁₁ is altered in a manner responsive to a transfer function of DUT 26 to become a forward DUT transmitted output signal _(O)Tx₁₂. Second probe 200 then extracts and terminates signal _(O)Tx₁₂. Similarly, any portion of signal _(I)Tx₁₁ that is reflected by or through DUT 26 returns to first probe 100 as forward DUT reflected output signal _(O)Tx₁₁. First probe 100 extracts and terminates signal _(O)Tx₁₁.

Port-2 transmission signal Tx_(2J), not being required, is never generated by transmitter 300. The path that would otherwise carry port-2 transmission signal Tx_(2J) is terminated within transmitter 300 and second probe 200.

Those skilled in the art will appreciate that a like process for reverse signal measurement (i.e., the signal to be measured passes from second probe 200, into second port 32, through DUT 26, out of first port 30, and into first probe 100) derives a reverse DUT input signal _(I)Tx₂₂ from port-2 transmission signal Tx₂₂, injects input signal _(I)Tx₂₂ into second port 32, alters input signal _(I)Tx₂₂ in a manner responsive to a reverse transfer function of DUT 26 to become a reverse DUT output signal _(O)Tx₂₁, and extracts and terminates output signal _(O)Tx₂₁. Port-1 transmission signal Tx_(1J) is, in this case, never generated by transmitter 300. The path that would otherwise carry port-1 transmission signal Tx_(1J) is terminated within transmitter 300 and first probe 100. Unless specifically indicated otherwise, forward signal measurement is assumed throughout this discussion.

Port-1 composite data and incident signals Dat₁₁ and Inc₁₁ are derived from forward DUT input signal _(I)Tx₁₁ within first probe 100, with Dat₁₁ being derived from _(I)Tx₁₁ through transmitted output _(O)Tx₁₁, propagated through high-frequency propagation medium 102 to become composite data and incident signals Dat₁₃ and Inc₁₃, and received by receiver 400 for processing. Similarly, a port-2 composite data signal Dat₂₂ is derived from forward DUT transmitted output signal _(O)Tx₁₂ within second probe 200, propagated through high-frequency propagation medium 202 to become a composite data signal Dat₂₃, and received by receiver 400 for processing.

FIGS. 2, 3, 4, and 5 show simplified schematic diagrams of transmitter 300, first probe 100, second probe 200, and receiver 400 and processor 500 of TES network analyzer 22 of FIG. 1, respectively, in accordance with a preferred embodiment of the present invention. The following discussion refers to FIGS. 1 through 5.

Port-1 composite data signal Dat₁₁ contains reflected output signal _(O)Tx₁₁, i.e., the reflected (reverse) component of forward DUT input signal _(I)Tx₁₁, at first port 30 of DUT 26 (i.e., in first probe 100 proximate first port 30 of DUT 26). Port-2 composite data signal Dat₂₂ contains the transmitted (forward) component of forward DUT transmitted output signal _(O)Tx₁₂ at second port 32. Similarly, port-1 composite incident signal Inc₁₁ contains the transmitted (forward) component of forward DUT input signal _(I)Tx₁₁ at first port 30. Signals Dat₁₁, Dat₂₂, and Inc₁₁ are composite signals. That is, each signal is formed by combining at least two signals at different frequencies.

Transmitter 300 develops several signals for transmission to first and second probes 100 and 200. Signal source 24 generates source signal _(S)Tx₃₃ at a given frequency and phase. Those skilled in the art will appreciate that source signal _(S)Tx₃₃ may be one of a variety of waveforms at any frequency. However, for exemplary purposes, source signal _(S)Tx₃₃ is assumed to be a swept-sinusoidal signal having, at any predetermined moment in time, a predetermined frequency f_(S) (not shown) having a predetermined phase Φ_(S) (not shown). Predetermined frequency f_(S) and phase Φ_(S) provide references against which other frequencies and phases within signal measurement system 20 are compared.

In the exemplary embodiment, signal source 24 provides a stable, high-frequency, swept CW signal _(S)Tx₃₃ to transmitter 300. For millimeter-wave applications, a relatively low frequency source may be used followed by frequency multiplication in the probe to avoid the excessive cable losses that occur at high frequencies. Signal source 24 is frequency locked to a system clock 502 and its frequency and output power are controlled by processor 500 via a system data and control bus 504 (e.g., a general-purpose information bus (GPIB)). The signal source may be, for example, a synthesizer which is synchronized (locked) to a 10 MHz clock. The 10 MHz clock also controls the triggering and sampling of the analog-to-digital converters (described hereinafter). When signal source 24 and receiver 400 are synchronized in this manner, the received signal will appear at substantially the same frequency as the transmitted signal, and there will be no leakage of the received signal into adjacent frequency measurement systems. This is particularly advantageous with the TES technique of the present invention.

Those skilled in the art will appreciate that it is also possible to apply the TES technique to systems with unsynchronized source signals. In such cases, it may be necessary to window the received signals in order to reduce leakage of the received frequency in adjacent frequencies. There may also be some loss in overall system accuracy.

Predetermined frequency f_(S) of signal source 24 is expressed as f _(S) =f _(IF) +f _(BF)+Δ_(LO),  (1) where f_(IF) denotes the frequency of an intermediate-frequency (IF) down-conversion source signal _(S)IF₃₃ (discussed hereinafter), f_(BF) denotes the frequency of a baseband down-conversion source signal _(S)BF₃₃ (discussed hereinafter), and Δ_(LO) denotes the approximate frequency for each of the local-oscillator signals (discussed hereinafter) associated with the frequency-offset mixers in probes 100 and 200. That is, f _(LPathP)=Δ_(LO)−δ_(LPathP).  (2)

In the above expression, f_(LPathP) denotes that we are referring to the frequency of a local-oscillator signal associated with one of the frequency-offset mixers in one of the signal-flow paths in the port-P probe 100 or 200. For example, f_(LPathP)=f_(LInc11) for a local-oscillator signal _(L)Inc₁₁ for the port-1 incident frequency-offset mixer (discussed hereinafter), while f_(LPathP)=f_(LComV22) for local-oscillator signal _(L)ComV₂₂ for the port-2 variable-common frequency-offset mixer (discussed hereinafter), etc. The parameter δ_(LPathP) denotes the difference frequency associated with one of the frequency-offset mixers. As is discussed hereinafter, the various offset signals are employed during signal normalization to remove the time-varying system errors.

Referring to FIGS. 2, 3, and 4, source signal _(S)Tx₃₃ is passed to transmitter 300. Source signal _(S)Tx₃₃ is level controlled by transmission signal attenuator 302 (FIG. 2) to produce a transmission signal Tx₃₃. Transmission signal Tx₃₃ is then passed to port-1 and port-2 transmission switches 304 and 306 (FIG. 2), respectively. Switches 304 and 306 work in conjunction with switches 106 (FIG. 3) and 206 (FIG. 4) and loads 308 and 310 (FIG. 2), 108 and 110 (FIG. 3), and 208 and 210 (FIG. 4) to determine the direction and loading of the signal through DUT 26.

For forward signal measurements, switches 304, 106, 306, and 206 are set to forward positions (the positions shown in FIGS. 2, 3, and 4). Transmission signal Tx₃₃ originates in transmitter 300 (FIG. 2). Transmission signal Tx₃₃ then passes through switch 304 to become port-1 transmission signal Tx₁₃. Signal Tx₁₃ is propagated over a port-1 transmission (Tx1) signal path 112 through high-frequency propagation medium 102 to become signal Tx₁₁ in first probe 100 (FIG. 3). Signal Tx₁₁ then passes through switch 106 to become forward DUT input signal _(I)Tx₁₁. Signal _(I)Tx₁₁ is propagated through DUT 26 to become forward DUT transmitted output signal _(O)Tx₁₂ in second probe 200 (FIG. 4) and is reflected by or through DUT 26 to become reflected output signal _(O)Tx₁₁ in first probe 100. In the exemplary embodiment, input and output signals _(I)Tx₁₁ and _(O)Tx₁₂ (i.e., the signal through DUT 26) are terminated by 50Ω loads 110 and 210 (since input signal _(I)Tx₁₁ is extracted from signal Tx₁₁ propagated over a 50Ω cable effectively in parallel with 50Ω load 110, input signal _(I)Tx₁₁ is loaded by 25Ω). Port-2 transmission signal Tx₂₃ is not generated, so switches 306 (FIG. 2) and 206 (FIG. 4) couple a port-2 transmission (Tx2) signal path 212 within high-frequency propagation medium 202 to dummy 50Ω loads 310 (FIG. 2) and 208 (FIG. 4). Those skilled in the art will appreciate that 50Ω impedances are exemplary only, and that other impedances may be used to meet the needs of specific circuits. The use of other impedances does not depart from the spirit of the present invention.

In second probe 200 (FIG. 4), 50Ω load 210 and switch 206, set to the forward position (as shown), provide a repeatable (static) load match for DUT 26, and isolate the time-varying effects of unused Tx2 signal path 212. If switch 206 were set to a reverse position (opposite the position shown), or if switch 206 and load 210 were not present, then the errors caused by a time-varying return loss in Tx2 signal path 212, carrying only a single signal, would not be removable using signal normalization techniques. For example, a portion of DUT output signal _(O)Tx₁₂ may be reflected by the time-varying impedance mismatch of Tx2 signal path 212. Part of this reflected signal may pass back through couplers 222 and 224 and DUT 26, where it may then be detected by a data coupler 124 (FIG. 3). Any time-varying cable return loss might thereby cause an error in the measured return loss of DUT 26. Such errors would be most severe when measuring low-return-loss DUTs.

Those skilled in the art will appreciate that other techniques may be used to remove these effects, including but not limited to the use of isolators, amplifiers, pads, etc. The use of other techniques does not depart from the spirit and intent of the present invention.

In first probe 100 (FIG. 3), switch 106 is set to the forward position (as shown) so that signal Tx₁₁ (i.e., signal _(I)Tx₁₁) may reach DUT 26. In the exemplary embodiment, Tx1 signal path 112 has an impedance of 50Ω. Since 50Ω load 110 is in parallel with Tx1 signal path 112, the net impedance is 25Ω. Switch 106 and 50Ω load 110 therefore contribute an insertion loss of 3.5 dB. It is possible to avoid this insertion loss if a wideband transformer is used to match the 50Ω impedance of Tx1 signal path 112.

Those skilled in the art will appreciate that switches 106 and 206 and loads 108, 110, 208, and 210 may be removed from Tx1 and Tx2 signal paths 112 and 212 if the time-varying load match is treated as an unknown parameter in the signal flow graph. The use of this and/or other techniques does not depart from the spirit of the present invention.

For reverse signal measurements, switches 304, 106, 306, and 206 are set to the reverse positions (opposite the positions shown in FIGS. 2, 3, and 4). Transmission signal Tx₃₃ (FIG. 2) passes through switch 306 to become port-2 transmission signal Tx₂₃. Signal Tx₂₃ is propagated over Tx2 signal path 212 through high-frequency propagation medium 202 to become signal Tx₂₂ in second probe 200 (FIG. 4). Signal Tx₂₂ then passes through switch 206 to become a reverse DUT input signal _(I)Tx₂₂ (not shown). Signal _(I)Tx₂₂ is propagated through DUT 26 in the reverse direction (i.e., from second port 32 to first port 30) to become a reverse DUT output signal _(O)Tx₂₁ (not shown) in first probe 100 (FIG. 3). Reverse input and output signals _(I)Tx₂₂ and _(O)Tx₂₁ (i.e., the reverse signal through DUT 26) are terminated by 50Ω loads 210 and 110. Port-1 transmission signal Tx₁₃ is not generated, so switches 304 (FIG. 2) and 106 (FIG. 3) couple Tx1 signal path 112 within high-frequency propagation medium 102 to dummy 50Ω loads 308 (FIG. 2) and 108 (FIG. 3). Again, unless specifically stated otherwise, forward signal measurement is presumed throughout this discussion.

Those skilled in the art will appreciate that the above transmission-signal scenarios are but a few of many techniques by which the measurement signal may be generated and propagated, and the use of other techniques does not depart from the spirit of the present invention.

In addition to transmission signals Tx₁₃ or Tx₂₃, transmitter 300 (FIG. 2) produces fixed-reference and variable-reference signals. A directional coupler 312 extracts an extracted source signal _(X)Tx₃₃ from source signal _(S)Tx₃₃. A first splitter 314 splits extracted source signal _(X)Tx₃₃ into port-1 and port-2 fixed-reference source signals _(S)RefF₁₃ and _(S)RefF₂₃, and a second splitter 316 splits extracted source signal _(X)Tx₃₃ into port-1 and port-2 variable-reference source signals _(S)RefV₁₃ and _(S)RefV₂₃, respectively. Reference source signals _(S)RefF₁₃, _(S)RefF₂₃, _(S)RefV₁₃, and _(S)RefV₂₃ are level controlled by variable attenuators 318, 320, 322, and 324, respectively, to become reference signals RefF₁₃, RefF₂₃, RefV₁₃, and RefV₂₃.

From transmitter 300 (FIG. 2), port-1 fixed-reference signal RefF₁₃ and variable-reference signal RefV₁₃ are propagated over fixed-reference (RefF1) and variable-reference (RefV1) signal paths 114 and 118 through high-frequency propagation medium 102 to become signals RefF₁₁ and RefV₁₁, respectively, in first probe 100 (FIG. 3). Then, a splitter 116 splits fixed-reference signal RefF₁₁ into link-common and fixed-common signals ComL₁₁ and ComF₁₁, and another splitter 120 splits variable-reference signal RefV₁₁ into variable-common and data-common signals ComV₁₁ and ComD₁₁. Similarly, from transmitter 300 (FIG. 2), port-2 fixed-reference signal RefF₂₃ and variable-reference signal RefV₂₃ are propagated over fixed-reference (RefF2) and variable-reference (RefV2) signal paths 214 and 218 through high-frequency propagation medium 202 to become signals RefF₂₂ and RefV₂₂, respectively, in second probe 200 (FIG. 4). Then, a splitter 216 splits fixed-reference signal RefF₂₂ into link-common and fixed-common signals ComL₂₂ and ComF₂₂, and another splitter 220 splits variable-reference signal RefV₂₂ into variable-common and data-common signals ComV₂₂ and ComD₂₂. Common signals ComL₁₁, ComF₁₁, ComV₁₁, and ComD₁₁ (FIG. 3), and ComL₂₂, ComF₂₂, ComV₂₂, and ComD₂₂ (FIG. 4) are signal-normalization signals.

For forward signal measurement, directional couplers 122 and 124 (FIG. 3) extract a raw incident signal _(R)Inc₁₁ and a reflected data signal _(R)Dat₁₁ from forward DUT input signal _(I)Tx₁₁ and reflected output signal _(O)Tx₁₁, respectively. Similarly, directional coupler 224 (FIG. 4) extracts a raw transmitted data signal _(R)Dat₂₂ from forward DUT output signal _(O)Tx₁₂. Directional coupler 222 (FIG. 4) and raw incident signal _(R)Inc₂₂ are not used in forward signal measurements.

In order to avoid potential cross-talk between the various signals and improve the isolation between the various signal-flow paths, mixing circuits (i.e., padded and filtered mixers) 126, 128, 130, 132, 134, and 136 (FIG. 3), and 228, 230, 232, 234, and 236 (FIG. 4) mix signals _(R)Inc₁₁, _(R)Dat₁₁, ComL₁₁, ComF₁₁, ComV₁₁, and ComD₁₁ (FIG. 3), and _(R)Dat₂₂, ComL₂₂, ComF₂₂, ComV₂₂, and ComD₂₂ (FIG. 4) at substantially predetermined frequency f_(S) (not shown) with stable, low-frequency signals _(L)Inc₁₁, _(L)Dat₁₁, _(L)ComL₁₁, _(L)ComF₁₁, _(L)ComV₁₁, and _(L)ComD₁₁ (FIG. 3), and _(L)Dat₂₂, _(L)ComL₂₂, _(L)ComF₂₂, _(L)ComV₂₂, and _(L)ComD₂₂ (FIG. 4) to produce offset signals _(O)Inc₁₁, _(O)Dat₁₁, _(O)ComL₁₁, (Link₁₁), _(O)ComF₁₁, _(O)ComV₁₁, and _(O)ComD₁₁ (FIG. 3), and _(O)Dat₂₂, _(O)ComL₂₂ (Link₂₂), _(O)ComF₂₂, _(O)ComV₂₂, and _(O)ComD₂₂ (FIG. 4), respectively, offset in frequency from predetermined frequency f_(S). That is, each of the mixing circuits produce an output signal whose constituent frequencies are substantially a sum and a difference of predetermined frequency f_(S) and the frequency of the particular low-frequency signal coupled to that mixing circuit. Filtration (discussed hereinafter) is used to select either the frequency-sum signal or the frequency-difference signal. In exemplary TES network analyzer 22, the frequency-difference signal is selected and the frequency-sum signal is attenuated. For example, the frequency of signal _(O)Inc₁₁ is the difference of the frequency of signal _(R)Inc₁₁ (substantially predetermined frequency f_(S)) less the frequency of local-oscillator signal _(L)Inc₁₁.

Referring to FIG. 2, signals _(L)PNC₃₃ (discussed hereinafter), _(L)Inc₁₃, _(L)Dat₁₃, _(L)ComL₁₃, _(L)ComF₁₃, _(L)ComV₁₃, _(L)ComD₁₃, _(L)PN₃₃ (discussed hereinafter), _(L)Inc₂₃, _(L)Dat₂₃, _(L)ComL₂₃, _(L)ComF₂₃, _(L)ComV₂₃, and _(L)ComD₂₃ are fixed, stable, low-frequency (e.g., approximately 100 kHz) signals from outputs 326, 328, 330, 332, 334, 336, 338, 340, 342, 344, 346, 348, 350, and 352 of a local-oscillator module 354. The frequency of each of the signals from local-oscillator module 354 is different (discussed in greater detail hereinafter) so that the frequencies of the resultant offset signals will be offset from predetermined frequency f_(S) by different amounts. Local-oscillator signals _(L)Inc₁₃, _(L)Dat₁₃, _(L)ComL₁₃, _(L)ComF₁₃, _(L)ComV₁₃, and _(L)ComD₁₃ (FIG. 2) are then propagated though low-frequency propagation medium 104 to first probe 100 (FIG. 3) to become signals _(L)Inc₁₁, _(L)Dat₁₁, _(L)ComL₁₁, _(L)ComF₁₁, _(L)ComV₁₁, and _(L)ComD₁₁, respectively. Similarly, local-oscillator signals _(L)Inc₂₃, _(L)Dat₂₃, _(L)ComL₂₃, _(L)ComF₂₃, _(L)ComV₂₃, and _(L)ComD₂₃ are propagated though low-frequency propagation medium 204 to second probe 200 (FIG. 4) to become signals _(L)Inc₂₂, _(L)Dat₂₂, _(L)ComL₂₂, _(L)ComF₂₂, _(L)ComV₂₂, and _(L)ComD₂₂, respectively.

The outputs of local-oscillator module 354 are also summed and monitored (discussed in greater detail hereinafter). Therefore, any phase drift associated with the outputs of local-oscillator module 354 is removed from the measurements. Thus, since transmission signals Tx₁₁, RefF₁₁, RefV₁₁, Tx₂₂, RefF₂₂, and RefV₂₂ are derived from source signal _(S)Tx₃₃, and reception signals _(O)Inc₁₁, _(O)Dat₁₁, _(O)ComL₁₁, _(O)ComF₁₁, _(O)ComV₁₁, and _(O)ComD₁₁, _(ODat) ₂₂, _(O)ComL₂₂, _(O)ComF₂₂, _(O)ComV₂₂, and _(O)ComD₂₂ are derived from the transmission signals, and since drift in source signal _(S)Tx₃₃ will be normalized out, the TES technique has no need of a precision, absolute-calibration source involving great expense and circuit complexity to minimize drift, but that may nevertheless still drift over time.

Port-1 offset link-common signal _(O)ComL₁₁ within first probe 100 (FIG. 2) is also port-1 link signal Link₁₁. Link signal Link₁₁ is combined into a composite link signal Link_(P1) (where P=1,2) by a combiner/splitter 138. Composite link signal Link_(P1) is propagated over a port-1 link (Link1) signal path 140 within high-frequency propagation medium 102 to become composite link signal Link_(P3) within receiver 400 (FIG. 5). Composite link signal Link_(P3) is propagated over a port-2 link (Link2) signal path 240 within high-frequency propagation medium 202 to become a composite link signal Link_(P2) within second probe 200 (FIG. 4). A port-1 link signal Link₁₂ is then split from composite link signal Link_(P2) by a combiner/splitter 238. Similarly, port-2 offset link-common signal _(O)ComL₂₂ within second probe 200 (FIG. 4) is also port-2 link signal Link₂₂. Link signal Link₂₂ is combined into composite link signal Link_(P2) by combiner/splitter 238. Composite link signal Link_(P2) is propagated over Link2 signal path 240 to become composite link signal Link_(P3) within receiver 400 (FIG. 5). Composite link signal Link_(P3) is propagated over Link1 signal path 140 to become composite link signal Link_(P1) within first probe 100 (FIG. 3). A port-2 link signal Link₂₁ is then split from composite link signal Link_(P1) by combiner/splitter 138. Thus, port-1 link signal Link₁₂ is present in second probe 200 and port-2 link signal Link₂₁ is present in first probe 100.

Within first probe 100 (FIG. 3), a combiner 142 combines offset incident signal _(O)Inc₁₁, offset variable-common signal _(O)ComV₁₁, offset fixed-common signal _(O)ComF₁₁, and port-2 link signal Link₂₁ (originally port-2 offset link-common signal _(O)ComL₂₂) to produce composite incident signal Inc₁₁. Composite incident signal Inc₁₁ is propagated over a port-1 incident (Inc1) signal path 144 within high-frequency propagation medium 102 to become composite incident signal Inc₁₃ in receiver 400 (FIG. 5). Additionally, a combiner 146 combines offset data signal _(O)Dat₁₁ and offset data-common signal _(O)ComD₁₁ to produce composite data signal Dat₁₁. Composite data signal Dat₁₁ is propagated over a port-1 data (Dat1) signal path 148 within high-frequency propagation medium 102 to become composite data signal Dat₁₃ in receiver 400.

Similarly, within second probe 200 (FIG. 4), a combiner 242 combines offset variable-common signal _(O)ComV₂₂, offset fixed-common signal _(O)ComF₂₂, and port-1 link signal Link₁₂ (originally port-1 offset link-common signal _(O)ComL₁₁) to produce a composite incident signal Inc₂₂. (Offset incident signal _(O)Inc₂₂ is not generated for forward signal measurements, and therefore not included in composite incident signal Inc₂₂.) Composite incident signal Inc₂₂ is propagated over a port-2 incident (Inc2) signal path 244 within high-frequency propagation medium 202 to become incident signal Inc₂₃ in receiver 400 (FIG. 5). Additionally, a combiner 246 combines offset data signal _(O)Dat₂₂ and offset data-common signal _(O)ComD₂₂ to produce composite data signal Dat₂₂. Composite data signal Dat₂₂ is propagated over a port-2 data (Dat2) signal path 248 within high-frequency propagation medium 202 to become composite data signal Dat₂₃ in receiver 400. Receiver 400 therefore receives composite signals Inc₁₃, Dat₁₃, Dat₂₃, and Inc₂₃ (signal Link_(P3) is a special case discussed herein)

Those skilled in the art will appreciate that the combination and propagation techniques described hereinbefore for composite signals Link_(P3), Inc₁₃, Dat₁₃, Dat₂₃, and Inc₂₃ are but some of the possible techniques, and that the use of other techniques does not depart from the spirit of the present invention. For example, multiple transmitter signals at distinct frequencies can be used to simultaneously measure the forward and reverse characteristics of DUT 26, thereby reducing the measurement time. In addition, a data path with minimal padding may be measured simultaneously with a data path with large padding for increased dynamic range.

In addition to the signal normalization discussed hereinbefore, the 12-term error model, well known to those skilled in the art, is used for static error suppression of first and second probes 100 and 200. After an initial calibration procedure utilizing the SOLT or other technique, the 12-term error model for probes 100 and 200 are computed. Once the 12-term error model has been computed, processor 500 may remove static systematic errors of probes 100 and 200 from future measurements. The combination of signal normalization and 12-term error modeling successfully suppresses repeatable errors from signal measurement system 20.

Whenever multiple signals reside on a single channel, noise can limit the overall accuracy with which measurements may be made. The effects of this noise are most severe at higher frequencies, or whenever small signals must simultaneously be measured in the presence of very large signals. Those skilled in the art will appreciate that the greatest preponderance of this noise is phase noise. Therefore, for the sake of simplicity, when this discussion refers to the suppression of phase noise, other noise will be assumed also to be suppressed.

TES network analyzer 22 provides phase-noise error suppression in two ways. First, the levels of the signals within each of composite signals Inc₁₃, Dat₁₃, Dat₂₃, and Inc₂₃ within receiver 400 (FIG. 5) are maintained within ±15 dB of one another. This is accomplished by including both fixed-reference signals RefF_(1J) and RefF_(2J) and variable-reference signals RefV_(1J) and RefV_(2J) into the design. In essence, as the level of composite data signals Dat₁₃ and Dat₂₃ change because of changes in the characteristics of DUT 26, a combination of gain ranging and variable common signal level adjustment is used to maintain the signal levels within the maximum ±15 dB range. This allows a measurement accuracy of 0.1 percent to be maintained for variation in DUT 26 over a 60 dB range, provided that the phase noise at the next-used frequency bin within processor 500 is less than −75 dBc (decibels relative to carrier).

Those skilled in the art will appreciate that the maintenance of the signals within each of composite signals Inc₁₃, Dat₁₃, Dat₂₃, and Inc₂₃ within ±15 dB is to some degree arbitrary. For example, TES network analyzer 22 could be designed to maintain the signals within each of composite signals Inc₁₃, Dat₁₃, Dat₂₃, and Inc₂₃ within ±20 dB of each other. While this would allow a measurement accuracy of 0.1 percent to be maintained for variations of DUT 26 over an 80 dB dynamic range, it would also lower the maximum phase noise requirement to −80 dBc. The use of other signal levels does not depart from the spirit of the present invention.

Referring to FIG. 5, receiver 400 passes signals received from first probe 100 (FIG. 3), second probe 200 (FIG. 4), and transmitter 300 (FIG. 2) through one unique and several substantially identical channels. These channels are Link channel 401 (the unique channel, discussed hereinafter), LOSum channel 402 (discussed hereinafter), Inc1 channel 403 and Dat1 channel 404 to process port-1 composite incident and data signals Inc₁₃ and Dat₁₃, respectively, Dat2 channel 405 and Inc2 channel 406 to process port-2 composite data and incident signals Dat₂₃ and Inc₂₃, respectively, and PN channel 407 (discussed hereinafter).

Gain ranging for channels 403, 404, 405, and 406 is effected through the use of variable-gain amplifiers 413, 414, 415, and 416 for composite signals Inc₁₃, Dat₁₃, Dat₂₃, and Inc₂₃, respectively. The TES technique is suited for gain ranging because a precise knowledge of the gain of each amplifier 413, 414, 415, and 416 is unnecessary since the amplifier's time-varying transfer function is removed via signal normalization within processor 500. Thus, the benefits of improved signal-to-noise ratios are provided by gain ranging, without the necessity of having to use precision components.

Whereas the first phase-noise error suppression technique discussed hereinbefore focuses on reducing the effects of phase noise by maintaining approximately equal levels among the signals within each of composite signals Inc₁₃, Dat₁₃, Dat₂₃, and Inc₂₃, the second phase-noise error suppression technique involves actually subtracting phase noise contamination from the measured signals. This technique is desirable at high frequencies where source signal _(S)Tx₃₃ (FIG. 2) does not meet the exemplary −75 dBc phase noise specification, or in applications where a better measurement accuracy or larger measurement dynamic range is desired.

The phase-noise-subtraction technique is based on the fact that the signals within each of composite signals Inc₁₃, Dat₁₃, Dat₂₃, and Inc₂₃ are derived from a common source signal _(S)Tx₃₃ (FIG. 2), thus exhibiting the same relative phase-noise spectra, and are all processed by receiver 400 substantially simultaneously in substantially identical channels 403, 404, 405, and 406, respectively. The technique requires that local-oscillator signals _(L)Inc₁₁, _(L)Dat₁₁, _(L)ComL₁₁, _(L)ComF₁₁, _(L)ComV₁₁, and _(L)ComD₁₁ (FIG. 3), and _(L)Inc₂₂, _(L)Dat₂₂, _(L)ComL₂₂, _(L)ComF₂₂, _(L)ComV₂₂, and _(L)ComD₂₂ (FIG. 4) contribute negligible phase noise and that the difference in time delay along each signal path is small relative to the measurement time.

A separate receiver channel, PN channel 407, is used to determine the phase-noise spectrum of source signal _(S)Tx₃₃ (FIG. 2). Processor 500 then subtracts the phase-noise spectrum of source signal _(S)Tx₃₃ from the phase-noise spectra of the signals within each of composite signals Inc₁₃, Dat₁₃, Dat₂₃, and Inc₂₂ in channels 403, 404, 405, and 406, respectively. A separate static calibration procedure (discussed hereinafter) is used to remove the frequency-dependent behavior from each of channels 403–407.

Within transmitter 300 (FIG. 2), a directional coupler 356 extracts an extracted source signal _(X)PN₃₃ from source signal _(S)Tx₃₃. A mixer 358 mixes signal _(X)PN₃₃ with a stable, low-frequency signal _(L)PN₃₃ from an output 340 of local-oscillator module 354 to produce a phase-noise reference signal PN₃₃. Transmitter 300 then passes signal PN₃₃ to receiver 400, where it is processed in PN channel 407. PN channel 407 is substantially identical to channels 403, 404, 405, and 406.

FIG. 2 depicts _(X)PN₃₃ extraction coupler 356 as a separate entity. This is done for clarity only. Those skilled in the art will appreciate that _(X)Tx₃₃ extraction coupler 312 serves substantially the same function (i.e., signals _(X)Tx₃₃ and _(X)PN₃₃ are substantially identical), and that combining couplers 312 and 356 into a single coupler does not depart from the spirit of the present invention.

FIG. 6 shows a simplified schematic of a generic mixing circuit 34 in accordance with a preferred embodiment of the present invention. The following discussion refers to FIG. 6.

Mixing circuit 34 (i.e., the mixer therein) produces an output signal whose frequency is either a sum or a difference of the frequencies of two input signals. Generic mixing circuit 34 has a first input signal In_(A) having a frequency f_(A), a second input signal In_(B) having a frequency f_(B), and an output signal Out having a frequency f_(O). Output frequency f_(O) is therefore one of a sum and a difference of input signal frequencies f_(A) and f_(B). That is, f_(O)=f_(A)+f_(B) or f_(O)=|f_(A)−f_(B)| (if high-order harmonics are ignored). If the difference frequency is chosen, output signal Out would be considered a down-conversion of input signal In_(A). This is typical of systems where f_(O) is a frequency much lower than f_(A). By being lower in frequency, output signal Out is more easily propagated though various frequency-variant propagation media with fewer errors and less distortion.

If, however, f_(B) is no greater than ten percent of f_(A), i.e., f_(B)≦0.1f_(A), then, while f_(O) does not equal f_(A), f_(O) is within ten percent (i.e., between 90 and 110 percent) of f_(A), i.e., 0.9f_(A)≦f_(O)≦1.1f_(A). Frequency f_(O) is close enough to frequency f_(A) so as to be affected in substantially the same manner by propagation through frequency-variant components and/or media of normal bandwidths. Carrying this argument further, if f_(B) is no greater than five percent of f_(A), i.e., f_(B)≦0.05f_(A), then f_(O) is between 95 and 105 percent of f_(A), i.e., 0.95f_(A)≦f_(O)≦1.05f_(A), and still further, if f_(B) is no greater than two percent of f_(A), i.e., f_(B)≦0.02f_(A), then f_(O) is between 98 and 102 percent of f_(A), i.e., 0.98f_(A)≦f_(O)≦1.02f_(A). It may readily be seen that as f_(B) becomes an ever-smaller percentage of f_(A), then f_(O) approaches f_(A), i.e., f_(O)≈f_(A). It is understood that f_(O) never equals f_(A), i.e., f_(O)≠f_(A).

Referring to FIGS. 2 through 6, to be offset rather than down-converted, the frequencies of local-oscillator signals _(L)Inc_(PJ), _(L)Dat_(PJ), _(L)ComL_(PJ), _(L)ComF_(PJ), _(L)ComV_(PJ), _(L)ComD_(PJ), and _(L)PN₃₃ (where P=1,2 and J=P,3) should be no greater than ten percent, preferably no greater than five percent, and more preferably no greater than two percent of predetermined frequency f_(S). In TES network analyzer 26 of the exemplary embodiment, source signal _(S)Tx₃₃ is a swept CW signal ranging in frequency from approximately 45 MHz to approximately 20 GHz, while the outputs of local-oscillator module 354 are as described in Table 1. These local-oscillator frequencies have been carefully chosen to avoid intermodulation distortion (IMD) between each of the signals produced by mixing each of the local-oscillator signals with the data, incident, reference, and link signals. A computer program was used to calculate all potential IMD products out to a specific order. Local-oscillator frequencies were rejected by the program if they produced an intermodulation component that potentially would interfere with the measured signals. Intermodulation components with 60 Hz and harmonics thereof were also calculated and rejected.

TABLE 1 Output Signal Freq. (Hz) 328 _(L)Inc₁₃ 99,610 330 _(L)Dat₁₃ 99,800 332 _(L)ComL₁₃ 99,320 334 _(L)ComF₁₃ 100,000 336 _(L)ComV₁₃ 99,450 338 _(L)ComD₁₃ 99,500 342 _(L)Inc₂₃ 99,330 344 _(L)Dat₂₃ 99,630 346 _(L)ComL₂₃ 99,360 348 _(L)ComF₂₃ 100,060 350 _(L)ComV₂₃ 99,720 352 _(L)ComD₂₃ 99,620

This produces a worst-case scenario where frequency f_(B) is 100.06 kHz (for local-oscillator signal _(L)ComF₂₃) and frequency f_(A) is 45 MHz (for input signal ComF₂₂), i.e., to the nearest 0.001 percent, f_(B)=0.00224f_(A). Therefore, frequency f_(O) (for output signal _(O)ComF₂₂) is only slightly offset from f_(A), i.e., 99.99776f_(A)≦f_(O)≦100.00224f_(A). Frequency f_(O) is affected in substantially the same manner as frequency f_(A) by propagation through frequency-variant components and/or media.

Those skilled in the art will appreciate that the frequencies specified in Table 1 for the local-oscillator signals are exemplary only, and that other frequencies may be used to suit the requirements of a given circuit. The use of such other frequencies does not depart from the spirit of the present invention.

Referring again to FIGS. 2 through 5, high-frequency propagation media 102 and 202 may affect signals Tx_(PJ), RefF_(PJ), and RefV_(PJ) during propagating from transmitter 300 to probes 100 and 200, respectively. Media 102 and 202 affect the simultaneous signals within each of composite signals Inc_(PJ) and Dat_(PJ) in substantially the same way during propagation from probes 100 and 200 to receiver 400. The signals within each of composite signals Link_(1J) and Link_(2J), propagating between probes 100 and 200, are similarly affected in a mutually identical manner. This equality of effect allows normalization to substantially suppress static and dynamic errors without recourse to high-precision devices and elaborate procedures.

In the exemplary embodiment, DUT 26 is coupled between first probe 100 and second probe 200, with first probe 100 connected to first port 30 and second probe 200 connected to second port 32. Conventionally, first port 30 and second port 32 are input and output ports, respectively, of DUT 26. In this way, when TES network analyzer 22 is set to measure a forward signal (i.e., a signal from first probe 100 to second probe 200 through DUT 26), the transmitted component (i.e., signal _(O)Tx₁₂) of signal _(I)Tx₁₁ passes through DUT 26 to second probe 200 in a normal input-to-output manner, with the reflected component (i.e., signal _(O)Tx₁₁) being reflected back to first probe 100.

Those skilled in the art will appreciate, however, that forward and reverse, and input and output, are relative terms, and that first and second probes 100 and 200 and first and second ports 30 and 32 may be interchanged at will without departing from the spirit of the present invention.

Transmitter 300 propagates transmission signals Tx_(1J), Ref_(1J), and RefV_(1J) to first probe 100, and transmission signals Tx_(2J), RefF_(2J), and RefV_(2J) to second probe 200 through high-frequency propagation media 102 and 202, respectively. First probe 100 propagates reception signals Dat_(1J) and Inc_(1J), and second probe 200 propagates reception signals Dat_(1J) and Inc_(1J) to receiver 400 through high-frequency propagation media 102 and 202, respectively. Also, first and second probes 100 and 200 reciprocally propagate signals Link_(1J) and Link_(2J) to second and first probes 200 and 100, respectively, through high-frequency propagation media 102 and 202. Transmitter 300 also propagates local-oscillator signals _(L)ComL_(1J), _(L)ComF_(1J), _(L)ComV_(1J), _(L)ComD_(1J), _(L)Inc_(1J), and _(L)Dat_(1J) to first probe 100 and local-oscillator signals _(L)ComL_(2J), _(L)ComF_(2J), _(L)ComV_(2J), _(L)ComD_(2J), _(L)Inc_(2J), and _(L)Dat_(2J) to second probe 200 through low-frequency propagation media 104 and 204, respectively.

High-frequency propagation media 102 and 202 may be bundles of flexible cables, with each cable providing one of signal paths 112, 114, 118, 140, 144, and 146 in high-frequency propagation medium 102 (FIGS. 2 and 3) and paths 212, 214, 218, 240, 244, and 246 in high-frequency propagation medium 202 (FIGS. 2 and 4). First probe 100 (FIG. 3) may be equipped with optional amplifiers 150, 152, 154, 156, and 158 on Tx1, RefF1, RefV1, Dat1, and Inc1 signal paths 112, 114, 118, 148, and 144, respectively. Similarly, second probe 200 (FIG. 4) may be equipped with optional amplifiers 250, 252, 254, 256, and 258 on Tx2, RefF2, RefV2, Dat2, and Inc2 signal paths 212, 214, 218, 248, and 244, respectively. Optional amplifiers 150, 152, 154, 156, 158, 250, 252, 254, 256, and 258 may be used when remote measurements of DUT 22 are required, and/or to model fiber optic or telemetry links. Dat1, Inc1, Dat2, and Inc2 optional amplifiers 156, 158, 256, and 258 may also be used to further isolate probes 100 and 200 from time-varying reflections associated with Dat1, Inc1, Dat2, and Inc2 signal paths 148, 144, 248, and 244, respectively.

Signal paths 112, 114, 118, 140, 148, and 144 (FIGS. 2 and 3), and 212, 214, 218, 240, 248, and 244 (FIGS. 2 and 4) in high-frequency propagation media 102 and 202 may be respectively characterized by time-varying scattering-matrices [S^(Tx1)], [S^(RefF1)], [S^(RefV1)], [S^(Link1)], [S^(Dat1)], [S^(Inc1)] (FIG. 3), and [S^(Tx2)], [S^(RefF2)], [S^(RefV2)], [S^(Link2)], [S^(Dat2)], and [S^(Inc2)] (FIG. 4). These scattering-matrices model the system drift effects that can occur along these paths, e.g., cable flexure errors, amplifier drift errors, and level variations resulting from changes in the settings of variable attenuators 302, 318, 320, 322, and/or 324 (FIG. 2), switches 304, 306 (FIG. 3), 106, and/or 206 (FIG. 4), and/or the gain of amplifiers 413, 414, 415, 416, and/or 417 (FIG. 5). Similar time-varying scattering-matrices are not necessary for low-frequency propagation media as cable flexure does not cause significant amplitude or phase changes at low frequencies. Furthermore, the mixer circuits are relatively insensitive to amplitude variations in the levels of the local-oscillator signals, thereby providing immunity to amplitude drift. Since the phases of the local-oscillator signals are monitored (discussed hereinafter), any phase drift in a local-oscillator signal is passed through the mixing circuit, and appears as a time-varying error in the mixed signal and is subtracted out in processor 500.

The circuitry in probes 100 and 200 serves two purposes. Directional couplers 122, 124, 222, and 224 are used to sample the incident, reflected, and transmitted signals associated with DUT 26. However, probes 100 and 200 also contain mixing circuits that allow the generation of signals _(O)Dat_(PP), _(O)Inc_(PP), _(O)ComL_(PP) (Link_(PP)), _(O)ComF_(PP), _(O)ComV_(PP), and _(O)ComD_(PP) whose frequencies are offset from predetermined frequency f_(S). That is, if higher-order harmonics are ignored, then the frequency-offset mixers produce signals at the following sum and difference frequencies: f _(S) +f _(LPathP) =f _(IF) +f _(BF)+2Δ_(LO)−δ_(LPathP); and  (3) f _(S) −f _(LPathP) =f _(IF) +f _(BF)+δ_(LPathP).  (4)

In the above expressions, f_(LPathP) denotes that we are referring to the frequency of a local-oscillator signal associated with one of the frequency-offset mixers in one of the signal-flow paths in the port-P probe 100 or 200 (where P=1,2). The parameter δ_(LPathP) denotes the difference frequency associated with one of the frequency-offset mixers. The technique of frequency offsetting signals _(O)Dat_(PP), _(O)Inc_(PP), _(O)ComL_(PP) (Link_(PP)), _(O)ComF_(PP), _(O)ComV_(PP), and _(O)ComD_(PP) eliminates signal crosstalk, improves the isolation between the various signal paths, and allows for the flow of multiple signals through common receiver circuitry.

The following discussion refers to FIGS. 2, 3, 4, and 6. Each of mixing circuits 126, 128, 130, 132, 134, 136, 226, 228, 230, 232, 234, 236, and 358 is equivalent to generic mixing circuit 34. That is, the discussion of generic mixing circuit 34 is simultaneously a discussion of each of mixing circuits 126, 128, 130, 132, 134, 136, 226, 228, 230, 232, 234, 236, and 358. Similarly, generic-circuit first input signal In_(A) is equivalent to any of mixing-circuit input signals _(R)Inc₁₁, _(R)Dat₁₁, ComL₁₁, ComF₁₁, ComV₁₁, ComD₁₁, _(R)Inc₂₂, _(R)Dat₂₂, ComL₂₂, ComF₂₂, ComV₂₂, ComD₂₂, and _(X)PN₃₃, generic-circuit second input signal In_(B) is equivalent to any of mixing-circuit local-oscillator signals _(L)Inc₁₁, _(L)Dat₁₁, _(L)ComL₁₁, _(L)ComF₁₁, _(L)ComV₁₁, _(L)ComD₁₁, _(L)Inc₂₂, _(L)Dat₂₂, _(L)ComL₂₂, _(L)ComF₂₂, _(L)ComV₂₂, _(L)ComD₂₂, and _(L)PN₃₃, and generic output signal Out is equivalent to any of mixing-circuit output (offset) signals _(O)Inc₁₁, _(O)Dat₁₁, _(O)ComL₁₁, _(O)ComF₁₁, _(O)ComV₁₁, _(O)ComD₁₁, _(O)Inc₂₂, _(O)Dat₂₂, _(O)ComL₂₂, _(O)ComF₂₂, _(O)ComV₂₂, _(O)ComD₂₂, and PN₃₃, respectively.

Generic mixing circuit 34 contains an input pad 36, an output pad 38, a low-frequency bandpass filter or diplexer 40, and a mixer 42. Diplexer 40 absorbs RF signals that may leak through on signal In_(B), thereby isolating mixer 42 from time-variations in signal In_(B), e.g., kinks in low-frequency propagation media 104 and 204, etc. Diplexer 40 also suppresses 60 Hz and harmonics thereof, and prevents any DC level in signal In_(B) from causing a time-varying attenuation of signal Out. It is also possible to make diplexer 40 a very narrow bandpass filter so that cross-talk between paths in low-frequency propagation media 104 and 204 may be eliminated, although crosstalk may also be virtually eliminated by using suitably shielded and/or twisted-pair cables in low-frequency propagation media 104 and 204. Since there is a large separation between the swept RF and local-oscillator frequencies, diplexers may also be employed on In_(A) and Out ports of generic mixing circuit 34 to filter out the low-frequency local-oscillator feedthrough.

Input and output pads 36 and 38 provide signal isolation, allow adjustment of signal levels, and provide impedance matching for mixer 42. This results in a reduction of the severity of intermodulation products generated by mixers 42.

The following discussion again refers to FIGS. 1 through 5. In order to make accurate measurements of the signal through DUT 26, it is necessary to establish accurate amplitude and phase references in probes 100 (FIG. 3) and 200 (FIG. 4). This is accomplished by producing local-probe link signals Link₁₁ and Link₂₂ in first and second probes 100 and 200, respectively. Signals Link₁₁ and Link₂₂ are then reciprocally propagated over a single bidirectional link (Link) signal path 140/240 (i.e., Link1 signal path 140 between first probe 100 and receiver 400, and Link2 signal path 240 between second probe 200 and receiver 400) through Link channel 401 (FIG. 5) to second and first probes 200 and 100 to become opposite-probe link signals Link₁₂ and Link₂₁, respectively. In opposite probes 200 and 100, opposite-probe link signals Link₁₂ and Link₂₁ are combined into composite incident signals Inc₂₂ and Inc₁₁, respectively, and passed to receiver 400.

Since link signals Link_(1J) and Link_(2J) (where J=1,2,3) are reciprocally propagated over single Link signal path 140/240, whatever propagation effect is imposed upon one link signal is also imposed upon the other link signal. Link signals Link_(1J) and Link_(2J) are used during signal normalization to measure the transmission coefficient through Link channel 401 (i.e., through Link signal path 140/240), thereby establishing a common tie between probes 100 and 200. Pads 160 and 162 (FIG. 3), and 260 and 262 (FIG. 4), as well as the output pads of mixing circuits 130 (FIG. 3) and 230 (FIG. 4), serve to isolate probes 100 and 200 from the effects of the time-varying return loss associated with Link signal path 140/240.

When DUT 26 is a reciprocal DUT (i.e., has the same transfer function for a forward and a reverse signal), some simplifications to the TES technique may be made. The path through DUT 26 is then also reciprocal, and may be used to establish the reference between probes 100 and 200. That is, the signal path through a reciprocal DUT 26 may also be used as a link channel (not shown). In TES network analyzer 22, the necessary probe-correlation information is carried on the incident and data signals Inc_(1J) and Dat_(1J) (FIG. 3), and Inc_(2J) and Dat_(2J) FIG. 4), and all the circuitry associated with fixed-reference signals RefF_(1J) (FIG. 3) and RefF_(2J) (FIG. 4) may then be removed (not shown).

Receiver 400 (FIG. 5) receives composite signals Inc₁₃, Dat₁₃, Dat₂₃, Inc₂₃, and PN₃₃ in channels 403, 404, 405, 406, and 407, respectively. Incident and data signals Inc₁₃, Dat₁₃, Dat₂₃, and Inc₂₃ are composite (multiple) signals, while phase-noise reference signal PN₃₃ is a sole (single) signal. In first probe 100 (FIG. 3), signals _(O)Inc₁₁, _(O)ComF₁₁, _(O)ComV₁₁, and Link₂₁ are combined into composite signal Inc₁₁ and propagated over Inc1 signal path 144 in high-frequency propagation medium 102 to become composite signal Inc₁₃ in receiver 400 (FIG. 5). Similarly, signals _(O)Dat₁₁, and _(O)ComD₁₁ are combined into composite signal Dat₁₁ and propagated over Dat1 signal path 148 in high-frequency propagation medium 102 to become composite signal Dat₁₃ in receiver 400. Likewise, in second probe 200 (FIG. 4), signals _(O)Inc₂₂, _(O)ComF₂₂, _(O)ComV₂₂, and Link₁₂ are combined into composite signal Inc₂₂ and propagated over Inc2 signal path 244 in high-frequency propagation medium 202 to become composite signal Inc₂₃ in receiver 400. Similarly, signals _(O)Dat₂₂, and _(O)ComD₂₂ are combined into composite signal Dat₂₂ and propagated over Dat2 signal path 248 in high-frequency propagation medium 202 to become composite signal Dat₂₃ in receiver 400.

In receiver 400 (FIG. 5), signals Inc₁₃, Dat₁₃, Dat₂₃, Inc₂₃, and PN₃₃ are gain ranged by variable-gain amplifiers 413, 414, 415, 416, and 417, respectively, as discussed hereinbefore (PN channel 407 processes signal PN₃₃ in exactly the same manner that channels 403, 404, 405, and 406 process composite signals Inc₁₃, Dat₁₃, Dat₂₃, and Inc₂₃). Gain-ranging amplifiers 413–417 also prevent the down-conversion local-oscillator feedthrough from reaching DUT 26, as may occur in conventional (i.e., prior-art) network analyzers.

Following gain ranging, gain-ranged signals Inc₁₃, Dat₁₃, Dat₂₃, Inc₂₃, and PN₃₃ pass though wideband filters (diplexers) 423, 424, 425, 426, and 427, respectively, before entering down-conversion mixers 433, 434, 435, 436, and 437. Diplexers 423, 424, 425, 426, and 427 suppress feed-through of signals _(L)Inc₁₁, _(L)ComF₁₁, _(L)ComV₁₁, or _(L)ComL₂₂ (FIG. 3) in composite signal Inc₁₃, signals _(L)Dat₁₁ or _(L)ComD₁₁ (FIG. 3) in composite signal Dat₁₃, signals _(L)Dat₂₂, or _(L)ComD₂₂ (FIG. 4) in composite signal Dat₂₃, signals _(L)Inc₂₂, _(L)ComF₂₂, _(L)ComV₂₂, or _(L)ComL₁₁ (FIG. 4) in composite signal Inc₂₃, and signal _(L)PN₃₃ (FIG. 2) in signal PN₃₃. Diplexers 423, 424, 425, 426, and 427 also provide good impedance matches for down-conversion mixers 433, 434, 435, 436, and 437, respectively.

Down-conversion mixers 433, 434, 435, 436, and 437 frequency multiply filtered signals Inc₁₃, Dat₁₃, Dat₂₃, Inc₂₃, and PN₃₃ with a swept-frequency, down-conversion source signal _(S)IF₃₃. Signal _(S)IF₃₃ of frequency f_(IF) is generated by an IF signal source (synthesizer) 438 driven by system clock 502 and synchronized with signal source 24 so as to maintain a constant difference in frequency between source signal _(S)Tx₃₃ and IF source signal _(S)IF₃₃. This difference is IF frequency f_(IF). The result is a down-conversion of signals Inc₁₃, Dat₁₃, Dat₂₃, Inc₂₃, and PN₃₃ into an IF incident signal _(I)Inc₁₃, an IF data signal _(I)Dat₁₃, an IF data signal _(I)Dat₂₃, an IF incident signal _(I)Inc₂₃, and an IF phase-noise reference signal _(I)PN₃₃, respectively. The frequencies of signals _(I)Inc₁₃, _(I)Dat₁₃, _(I)Dat₂₃, _(I)Inc₂₃, and _(I)PN₃₃ are offset from IF frequency f_(S)±f_(IF) to the same degree signals Inc₁₃, Dat₁₃, Dat₂₃, Inc₂₃, and PN₃₃ are offset from predetermined frequency f_(S). This produces upper-sideband signals at f _(S) +f _(LPathP) +f _(IF)=2f _(IF) +f _(BF)+2Δ_(LO)−δ_(LPathP), and  (5) f _(S) −f _(LPathP) +f _(IF)=2f _(IF) +f _(BF)+δ_(LPathP),  (6) and lower sidebands at f _(S) +f _(LPathP) −f _(IF) =f _(BF)+2δ_(LO)−δ_(LPathP), and   (7) f _(S) −f _(LPathP) −f _(IF) =f _(BF)+δ_(LPathP).  (8)

Those skilled in the art will appreciate that while TES network analyzer 22 employs fundamental mixing for the IF down-conversion, other techniques, such as harmonic mixing or sampling, may also be employed without departing from the spirit of the present invention. Furthermore, in millimeter-wave applications, where cable losses are prohibitively large, it may be desirable to place IF down-conversion mixers 433–436 in probes 100 and 200.

Signals _(I)Inc₁₃, _(I)Dat₁₃, _(I)Dat₂₃, _(I)Inc₂₃, and _(I)PN₃₃ are then passed through pads 443, 444, 445, 446, and 447, very-narrow-band crystal filters 453, 454, 455, 456, and 457, and down-conversion mixers 463, 464, 465, 466, and 467.

Pads 443–447 isolate mixers 433–437 from reflections resulting from impedance mismatches associated with crystal filters 453–457.

Very-narrow-band crystal filters 453–457 pass the difference signals from mixing circuits 126, 128, 130, 132, 134, 136, 226, 228, 230, 232, 234, 236, and 358. Since the frequencies of local-oscillator signals _(L)Inc_(PJ), _(L)Dat_(PJ), _(L)ComL_(PJ), _(L)ComF_(PJ), _(L)ComD_(PJ), and _(L)PN₃₃ are approximately 100 kHz, mixing circuits 126, 128, 130, 132, 134, 136, 226, 228, 230, 232, 234, 236, and 358 produce dual output signals whose frequencies are approximately 100 kHz above and below predetermined frequency f_(S) (not shown). After down-conversion by mixers 463–467, those frequencies are approximately 100 kHz above and below intermediate frequency f_(S)−f_(IF) (not shown). Filters 453–457 remove those signals whose frequencies are defined in Equations (5)–(7), leaving only those signals whose frequencies are given in Equation (8).

The outputs of crystal filters 453–457 flow into the inputs of down-conversion mixers 463–467. Mixers 463–467 are driven by fixed-frequency baseband source signal _(S)BF₃₃ of frequency f_(BF) derived from a baseband signal source (divider) 468 driven by system clock 502. In this manner, IF signals _(I)Inc₁₃, _(I)Dat₁₃, _(I)Dat₂₃, _(I)Inc₂₃, and _(I)PN₃₃ are converted into baseband signals _(B)Inc₁₃, _(B)Dat₁₃, _(B)Dat₂₃, _(B)Inc₂₃, and _(B)PN₃₃. The sum and difference frequencies for mixers 463–467 exist at f _(S) −f _(LPathP) −f _(IF) +f _(BF)=2f _(BF)+δ_(LPathP), and  (9) f _(S) −f _(LPathP) −f _(IF) −f _(BF)=δ_(LPathP).  (10)

In addition to baseband signals _(B)Inc₁₃, _(B)Dat₁₃, _(B)Dat₂₃, _(B)Inc₂₃, and _(B)PN₃₃, receiver 400 (FIG. 5) also receives a local-oscillator summing signal _(L)Sum₃₃ from transmitter 300 (FIG. 2). Within local-oscillator module 354 (FIG. 2), local-oscillator outputs 328, 330, 332, 334, 336, 338, 342, 344, 346, 348, 350, and 352, i.e., all the individual local-oscillator signals _(L)ComL₁₃, _(L)ComF₁₃, _(L)ComV₁₃, _(L)ComD₁₃, _(L)Inc₁₃, _(L)Dat₁₃, _(L)ComL₂₃, _(L)ComF₂₃, _(L)ComV₂₃, _(L)ComD₂₃, _(L)Inc₂₃, and _(L)Dat₂₃, respectively, are summed together in a summer 360 to produce summing signal _(L)Sum₃₃. Summing signal _(L)Sum₃₃ is processed by receiver 400 in LOSum channel 402.

Local-oscillator module 354 (FIG. 2) generates local-oscillator output signals _(L)PNC₃₃, _(L)ComL₁₃, _(L)ComF₁₃, _(L)ComV₁₃, _(L)ComD₁₃, _(L)Inc₁₃, _(L)Dat₁₃, _(L)PN₃₃, _(L)ComL₂₃, _(L)ComF₂₃, _(L)ComV₂₃, _(L)ComD₂₃, _(L)Inc₂₃, and _(L)Dat₂₃ for TES network analyzer 22. Since the frequencies for these output signals are relatively low (on the order of 100 kHz), local-oscillator module 354 may utilize inexpensive dividers or digital-to-analog converters to generate the output signals from system clock 502 under control of processor 500. Summer 360 sums all probe output signals together to produce local-oscillator summing signal _(L)Sum₃₃. Summer 360 may consist of, for example, 10 kΩ resistors from each of outputs 328, 330, 332, 334, 336, 338, 342, 344, 346, 348, 350, and 352 to a single connection point. The 10 kΩ resistors will not measurably alter the 50Ω impedance of the _(L)ComL₁₃, _(L)ComF₁₃, _(L)ComV₁₃, _(L)ComD₁₃, _(L)Inc₁₃, _(L)Dat₁₃, _(L)ComL₂₃, _(L)ComF₂₃, _(L)ComV₂₃, _(L)ComD₂₃, _(L)Inc₂₃, and _(L)Dat₂₃ signal paths, and will provide a small, but easily measurable, sample of the local-oscillator output signals. Since all of the local-oscillator output signals for probes 100 and 200 are summed together with this circuit, and since we only need relative phases of the local-oscillator output signals, any subsequent phase shifts of these combined signals will have no effect.

Those skilled in the art will appreciate that the above techniques for generating local-oscillator signals _(L)PNC₃₃, _(L)ComL₁₃, _(L)ComF₁₃, _(L)ComV₁₃, _(L)ComD₁₃, _(L)Inc₁₃, _(L)Dat₁₃, _(L)PN₃₃, _(L)ComL₂₃, _(L)ComF₂₃, _(L)ComV₂₃, _(L)ComD₂₃, _(L)Inc₂₃, and _(L)Dat₂₃, and for summing the probe signals into signal _(L)Sum₃₃, are exemplary. The use of other techniques does not depart from the spirit of the present invention.

Referring again to FIG. 5, signals _(L)Sum₃₃, _(B)Inc₁₃, _(B)Dat₁₃, _(B)Dat₂₃, _(B)Inc₂₃, and _(B)PN₃₃ pass through amplifiers 472, 473, 474, 475, 476, and 477 and low-pass filters 482, 483, 484, 485, 486, and 487 before being fed into high-resolution analog-to-digital converters (ADCs) 492, 493, 494, 495, 496, and 497.

Amplifiers 472, 473, 474, 475, 476, and 477 may drift slightly with temperature and/or over time. Any such drift may cause a slight shift in phase of signals _(L)Sum₃₃, _(B)Inc₃₃, _(B)Dat₁₃, _(B)Dat₂₃, _(B)Inc₂₃, and _(B)PN₃₃ passing through amplifiers 472, 473, 474, 475, 476, and 477, respectively. However, signal _(L)Sum₃₃, being a sum of local-oscillator signals _(L)ComL₁₃, _(L)ComF₁₃, _(L)ComV₁₃, _(L)ComD₁₃, _(L)Inc₁₃, _(L)Dat₁₃, _(L)ComL₂₃, _(L)ComF₂₃, _(L)ComV₂₃, _(L)ComD₂₃, _(L)Inc₂₃, and _(L)Dat₂₃, is digitized separately and later used to correct the measured data for any relative phase drift between the local-oscillator sources. Drift in the other signals _(B)Inc₁₃, _(B)Dat₁₃, _(B)Dat₂₃, _(B)Inc₂₃, and _(B)PN₃₃ is removed via signal normalization.

Low-pass filters 482–487 are designed to pass the desired signals while providing anti-alias protection for ADCs 492–497. ADCs 492–497 are used to digitally sample the relatively low frequency waveforms of signals _(L)Sum₃₃, _(B)Inc₁₃, _(B)Dat₁₃, _(B)Dat₂₃, _(B)Inc₂₃, and _(B)PN₃₃, convert those waveforms from the analog to the digital domain, and produce digital signals _(D)Sum₃₃, _(D)Inc₁₃, _(D)Dat₁₃, _(D)Dat₂₃, _(D)Inc₂₃, and _(D)PN₃₃. Signals _(D)Sum₃₃, _(D)Inc₁₃, _(D)Dat₁₃, _(D)Dat₂₃, _(D)Inc₂₃, and _(D)PN₃₃, being digital, are logical signals which pass over logical paths in channels 402–407 to processor 500. These logical paths are represented at the bottom of receiver 400 (FIG. 5) and in processor 500, (FIG. 7, discussed hereinafter) as dotted lines. Physically, signals _(D)Sum₃₃, _(D)Inc₁₃, _(D)Dat₁₃, _(D)Dat₂₃, _(D)Inc₂₃, and _(D)PN₃₃ are passed to processor 500 over system data and control bus 504.

Those skilled in the art will appreciate that the scenario discussed hereinbefore for receiver 400 is but one of a plurality of such scenarios, and that the use of another scenario does not depart from the spirit of the present invention. For example, since all the individual signals making up the composite signals Inc₁₃, Dat₁₃, Dat₂₃, and Inc₂₃ exist at distinct frequencies, the signals associated with two or more of receiver channels 403–406 may be combined together in a single receiver channel to reduce receiver costs.

The following discussion refers to FIGS. 2 and 5. The basic idea behind the phase-noise subtraction technique mentioned hereinbefore is to determine the level of the phase noise relative to source signal _(S)Tx₃₃ in dedicated PN channel 407 containing single (sole) signal PN₃₃. This information is then used by processor 500 to extract the phase-noise contamination in channels 403, 404, 405, and 406 from each of the signals within composite signals Inc₁₃, Dat₁₃, Dat₂₃, and Inc₂₃, respectively. It is assumed that the differences between the time delays associated with signal propagation along channels 403–407 are small in comparison to the data acquisition time. Furthermore, it is assumed that all of the signals are converted substantially simultaneously by ADCs 493–497. Under these assumptions, the phase-noise subtraction technique works since signals Inc₁₃, Dat₁₃, Dat₂₃, Inc₂₃, and PN₃₃ are all derived from source signal _(S)Tx₃₃, i.e., from high-frequency synthesized signal source 24, which is the cause of the phase noise contamination. IF and baseband down-conversion sources 438 and 468 may also contribute appreciable amounts of phase noise to the system. However, this phase noise contamination may be monitored and removed since it is common to signals Inc₁₃, Dat₁₃, Dat₂₃, Inc₂₃, and PN₃₃, i.e., IF and baseband source signals _(S)IF₃₃ and _(S)BF₃₃ are each applied to channels 403–407. The relatively low-frequency offset local-oscillator sources 328, 330, 332, 334, 336, 338, 340, 342, 344, 346, 348, 350, and 352, which are not common to all the measured signals, contribute negligible levels of phase-noise contamination, and may therefore be ignored.

An initial calibration procedure is used to remove static (systematic) errors associated with receiver channels 403–407 for the phase-noise subtraction technique discussed hereinbefore. Within transmitter 300 (FIG. 2), a mixer 362 and directional couplers 364 and 366 are used to inject moderate sized, simulated “phase-noise” sidebands into system 20. Directional coupler 364 extracts a signal _(X)PNC₃₃ from source signal _(S)Tx₃₃. Mixer 362 then mixes local-oscillator phase-noise-calibration signal _(L)PNC₃₃ with extracted signal _(X)PNC₃₃ to produce an injection signal _(I)PNC₃₃. A low-frequency bandpass filter 368 is used to remove RF components from local-oscillator signal _(L)PNC₃₃. Injection signal _(I)PNC₃₃ consists substantially of a sum and a difference of predetermined frequency f_(S) and the frequency of local-oscillator signal _(L)PNC₃₃. That is, injection signal _(I)PNC₃₃ consists of two “sidebands” of predetermined frequency f_(S). Directional coupler 366 then injects injection signal _(I)PNC₃₃ into source signal _(S)Tx₃₃ to produce a phase-noise source signal _(S)PN₃₃, hence producing the calibration sidebands in all signals in system 20. Phase-noise source signal _(S)PN₃₃ is substantially source signal _(S)Tx₃₃ plus the upper and lower sidebands from injection signal _(I)PNC₃₃, which sidebands simulate quantifiable phase noise with which to calibrate PN channel 407 against the sidebands in a signal in one of the other channels. Local-oscillator calibration signal _(L)PNC₃₃ may be varied to produce an optimal set of sidebands for calibration of any signal within system 20. Variable attenuators 302, 318, 320, 322, and 324 may be used to “turn-off” all the other signals within system 20 except the signal being calibrated during each step of the phase-noise subtraction technique calibration procedure.

Those skilled in the art will appreciate that the above phase-noise calibration technique is exemplary, and that the use of another technique does not depart from the spirit of the present invention.

FIG. 7 shows a simplified schematic diagram of processor 500 of TES network analyzer 22 in accordance with a preferred embodiment of the present invention. The following discussion refers to FIG. 7.

Processor 500 controls the signal routing, collection of the data, analysis of the data, and display of the processed data. System data and control bus 504 serves as the interface between processor 500 and other circuits of analysis unit 28. Additionally, processor 500 may generate system clock 502 in a system clock source 506. System clock 502 is then distributed to other portions of signal measurement system 20.

Those skilled in the art will appreciate that not all connections of system clock 502 and system data and control bus 504 are depicted in the Figures, and that those connections that are depicted are exemplary only. It will be further appreciated that there is no requirement that system clock source 506 be a component of processor 500, and that the use of other techniques of producing system clock 502 do not depart from the spirit of the present invention.

Processor 500 controls several functions of signal measurement system 20. For example, processor 500 sends commands to swept-frequency signal source 24 and IF signal source 438 that control the type of frequency sweep, the range of the frequency sweep, and the level of source signal _(S)Tx₃₃, processor 500 sets switches 304 and 306 (FIG. 2), 106 (FIG. 3), and 206 (FIG. 4) to select between forward DUT input signal _(I)Tx₁₁ and reverse DUT input signal _(I)Tx₂₂ as the signal to be measured, processor 500 sets variable attenuators 302, 318, 320, 322, and 324 (FIG. 2) to establish the levels of signals Tx₃₃, RefF₁₃, RefF₂₃, RefV₁₃, and RefV₂₃, respectively, processor 500 controls local-oscillator module 354 (FIG. 2) so as to establish the frequencies of signals _(L)PNC₃₃, _(L)ComL₁₃, _(L)ComF₁₃, _(L)ComV₁₃, _(L)ComD₁₃, _(L)Inc₁₃, _(L)Dat₁₃, _(L)PN₃₃, _(L)ComL₂₃, _(L)ComF₂₃, _(L)ComV₂₃, _(L)ComD₁₃, _(L)Inc₂₃, and _(L)Dat₂₃, processor 500 sets the gain of variable-gain amplifiers 413, 414, 415, and 416 (FIG. 5) so that the signals within each of composite signals Inc₁₃, Dat₁₃, Dat₂₃, Inc₂₃ are within ±15 db of each other, and processor 500 also establishes the sample time, sampling rate, and data-collection starting time of ADCs 492–493 (FIG. 5).

Once ADCs 492–497 (FIG. 5) have produced the digitized signals _(D)LSum₃₃, _(D)Inc₁₃, _(D)Dat₁₃, _(D)Dat₂₃, _(D)Inc₂₃, and _(D)PN₃₃, then these signals are sent through the system control/data bus 504 to processor 500. For purposes of clarity, FIG. 7 depicts signals _(D)LSum₃₃, _(D)Inc₁₃, _(D)Dat₁₃, _(D)Dat₂₃, _(D)Inc₂₃, and _(D)PN₃₃ as continuing to propagate along channels LOSum 402, Inc1 403, Dat1 404, Dat2 405, Inc2 406, and PN 407. Those skilled in the art will appreciate that this is a convention of the figures only, and that all signal propagate through bus 504 to processor 500, then propagate within processor 500 in a manner well-known for processors.

Once in processor 500, a conventional signal-averaging or decimating routine 508 is used to preprocess the signals. After preprocessing, a Fast Fourier Transform (FFT) routine 510 is used to transform signals _(D)LSum₃₃, _(D)Inc₁₃, _(D)Dat₁₃, _(D)Dat₂₃, _(D)Inc₂₃, and _(D)PN₃₃ from the time domain to the frequency domain and distribute the resultant _(T)LSum₃₃, _(T)Inc₁₃, _(T)Dat₁₃, _(T)Dat₂₃, _(T)Inc₂₃, and _(T)PN₃₃ signals into a plurality of frequency bins. A normalization routine 512 then extracts the desired, complex, frequency-domain signals _(T)LSum₃₃, _(T)Inc₁₃, _(T)Dat₁₃, _(T)Dat₂₃, _(T)Inc₂₃, and _(T)PN₃₃ from their respective frequency bins. A compensation routine 514 then uses normalized signals _(N)LSum₃₃, _(N)Inc₁₃, _(N)Dat₁₃, _(N)Dat₂₃, _(N)Inc₂₃, and _(N)PN₃₃ with the hereinbefore described static, dynamic, and phase-noise error suppression techniques to produce the desired error corrected S-parameters for DUT 26.

Once the S-parameters have been computed, processor 500 may display these data to the user (not shown) or otherwise effect desired actions (not shown). For exemplary TES network analyzer 22, the user has control over the type of display (i.e., log magnitude, linear magnitude, real and imaginary parts, or Smith chart representations), the frequency span, the number of frequency points within the frequency span, etc.

Since the digitized, time-domain waveforms are processed using FFT routine 510, all of the sources and samplers must be precisely synchronized. This synchronization is achieved using system clock 502. System clock 502 is used to synchronize signal source (swept-frequency synthesizer) 24 (FIG. 1), frequency-offset local-oscillator outputs (digital-to-analog converters) 326, 328, 330, 332, 334, 336, 338, 340, 342, 344, 346, 348, 350, and 352 (FIG. 2), IF down-conversion signal source 438 (FIG. 5), baseband down-conversion signal source 468, and ADCs (high-resolution analog-to-digital converters) 492, 493, 494, 495, 496, and 497 (FIG. 5).

A better understanding of TES network analyzer 22 and the signal normalization techniques of the TES technique can be gained by studying the signal-flow paths. The time-varying amplitude and phase errors associated with cable flexure and temperature drift are modeled by scattering-matrices (S-matrices) [S^(Tx1)], [S^(RefF1)], [S^(RefV1)], [S^(Link1)], [S^(Dat1)], and [S^(Inc1)] for port-1 high-frequency propagation medium 102 (FIG. 3), and [S^(Tx2)], [S^(RefF2)], [S^(RefV2)], [S^(Link2)], [S^(Dat2)], and [S^(Inc2)] for port-2 high-frequency propagation medium 202 (FIG. 4). These matrices may be generalized as [S^(PathP)], where the superscript Path indicates the signal-flow path of the specific transmitter signal Tx_(1J) or TX_(2J), receiver signal Dat_(1J), Dat_(2J), Inc_(1J), or Inc_(2J), reference signal RefF_(1J), RefF_(2J), RefV_(1J), or RefV_(2J), or link signal Link_(1J) or Link_(2J), where P denotes the port (P=1,2), and the square brackets indicate that we are referring to a matrix. Since these paths either contain a single signal that is later split into data and normalization signals with different frequencies, or they contain two or more simultaneous signals at different frequencies, signal normalization techniques can be used to remove both the static and dynamic errors associated with components in these paths from the measured S-parameters.

When discussing the individual S-parameters associated with DUT 26, the square brackets are dropped and subscripts are added to denote the transmitter and receiver ports (Table 2). Note that the time-varying characteristics of the optional transmission and reception amplifiers 150, 152, 154, 156, and 158 (FIG. 3), and 250, 252, 254, 256, and 258 (FIG. 4), and other signal-conditioning circuitry can also be included in these time-varying S-matrices. Optional Dat1, Inc1, Dat2, and Inc2 amplifiers 156 and 158, and 256 and 258 can be used to isolate probes 100 (FIG. 3) and 200 (FIG. 4) from reflections in Dat1 and Inc1 signal paths 148 and 144, and Dat2 and Inc2 signal paths 248 and 244, respectively. Optional transmission and reception amplifiers 150, 152, 154, 156, 158, 250, 252, 254, 256, and 258 may also be desirable in applications involving measurements of remote or high-power DUTs. In some applications, it may even be desirable to employ fiber optic or telemetry links to connect between transmitter 300 and/or receiver 400 and calibration probes 100 and 200. These effects can also be included in the time-varying S-matrices.

FIG. 8 shows a signal-flow graph of the signal-flow path of composite signal Dat₁₁ in accordance with a preferred embodiment of the present invention. The following discussion refers to FIGS. 2, 3, 4, and 8.

Composite signals Inc_(1J), Dat_(1J), Dat_(2J), and Inc_(2J) are made up of independent measured signals. These measured signals maintain their independence as combined signals Inc_(1J), Dat_(1J), Dat_(2J), and Inc_(2J) are propagated through and processed by receiver 400 (FIG. 5). For port-1 signals, combined signal Inc_(1J) contains independent measured signals _(M)Inc₁₃, _(M)ComV₁₃, _(M)ComF₁₃, and _(M)ComL₂₃ (_(M)Link₂₃) (none shown) derived from signals _(O)Inc₁₁, _(O)ComV₁₁, _(O)ComF₁₁, and _(O)ComL₂₁ (Link₂₁) combined by combiner 142, and combined signal Dat_(1J) contains independent measured signals _(M)Dat₁₃ and _(M)ComD₁₃ (none shown, save _(M)Dat₁₃ is an exemplary signal in FIG. 8), derived from signals _(O)Dat₁₁ and _(O)ComD₁₁ combined by combiner 146 in first probe 100 (FIG. 3). Similarly, for port-2 signals, combined signal Inc_(2J) contains independent measured signals _(M)Inc₂₃, _(M)ComV₂₃, _(M)ComF₂₃, and _(M)ComL₁₃ (_(M)Link₁₃) (none shown) derived from signals _(O)Inc₂₂, _(O)ComV₂₂, _(O)ComF₂₂, and _(O)ComL₁₂ (Link₁₂) combined by combiner 242, and combined signal Dat_(2J) contains independent measured signals _(M)Dat₂₃ and _(M)ComD₂₃ (none shown) derived from signals _(O)Dat₂₂ and _(O)ComD₂₂ combined by combiner 246 in second probe 200 (FIG. 4).

Assuming switches 304 and 306 (FIG. 2), 106 (FIG. 3), and 206 (FIG. 4) are set for forward DUT signal measurements, i.e., as shown in FIGS. 2, 3, and 4, then the signal-flow graph in FIG. 8 can be used to better understand signal flow along a signal-flow path of signal Dat_(1J). Other signal-flow paths in TES network analyzer 22 are similar. Thus, once the signal-flow analysis has been carried out for the path of signal Dat_(1J), then these results can be directly extended to the paths of other signals by replacing the variables listed in the _(M)Dat₁₃ row in Table 2 with the corresponding variables for another flow path from another row in Table 2.

TABLE 2 FFT Signal Transmitter Circuitry Probe Circuitry and DUT Receiver Circuitry _(M)Dat₁₃ Γ_(S) ^(Tx1) [S^(Tx1)] Γ_(L) ^(Tx1) {tilde over (S)}_(11M)P_(ODat11)exp(−j∠A_(LDat11)) Γ_(S) ^(Odat 1) [S^(Odat1)] Γ_(L) ^(Odat1) _(M)ComD₁₂ Γ_(S) ^(RefV 1) [S^(RefV1)] Γ_(L) ^(RefV 1) P_(OComD11)exp(−j∠A_(LComD11)) Γ_(S) ^(Odat 1) [S^(Odat1)] Γ_(L) ^(Odat1) _(M)Inc₁₃ Γ_(S) ^(Tx1) [S^(Tx1)] Γ_(L) ^(Tx1) P_(OInc11) exp(−j∠A_(LInc11)) Γ_(S) ^(OInc1) [S^(OInc1)] Γ_(L) ^(OInc1) _(M)ComV₁₃ Γ_(S) ^(RefV 1) [S^(RefV1)] Γ_(L) ^(RefV 1) P_(OComV11) exp(−j∠A_(LComV11)) Γ_(S) ^(OInc1) [S^(OInc1)] Γ_(L) ^(OInc1) _(O)ComF₁₃ Γ_(S) ^(RefF 1) [S^(RefF1)] Γ_(L) ^(RefF1) P_(OComF11) exp(−j∠A_(LComF11)) Γ_(S) ^(OInc1) [S^(OInc1)] Γ_(L) ^(OInc1) _(M)Link₂₃ Γ_(S) ^(RefF 2) [S^(RefF2)] Γ_(L) ^(RefF 2) T^(Link) P_(Link22) exp(−j∠A_(LComL22)) Γ_(S) ^(OInc1) [S^(OInc1)] Γ_(L) ^(OInc1) _(M)Dat₂₃ Γ_(S) ¹ [S^(Tx1)] Γ_(L) ^(Tx1) {tilde over (S)}_(21M)P_(ODat22) exp(−j∠A_(LDat22)) Γ_(S) ^(Odat2) [S^(Odat2)] Γ_(L) ^(Odat2) _(M)ComD₂₃ Γ_(S) ^(RefV 2) [S^(RefV2)] Γ_(L) ^(RefV 2) P_(OComD22) exp(−j∠A_(LComD22)) Γ_(S) ^(Odat2) [S^(Odat2)] Γ_(L) ^(Odat2) _(M)Inc₂₃ NA NA NA NA NA NA NA _(M)ComV₂₃ Γ_(S) ^(RefV 2) [S^(RefV2)] Γ_(L) ^(RefV 2) P_(OComV22) exp(−j∠A_(LComV22)) Γ_(S) ^(OInc2) [S^(OInc2)] Γ_(L) ^(OInc2) _(M)ComF₂₃ Γ_(S) ^(RefF 2) [S^(RefF2)] Γ_(L) ^(RefF 2) P_(OComF22) exp(−j∠A_(LComF22)) Γ_(S) ^(OInc2) [S^(OInc2)] Γ_(L) ^(OInc2) _(M)Link₁₃ Γ_(S) ^(RefF 1) [S^(RefF1)] Γ_(L) ^(RefF 1) T^(Link) P_(Link11) exp(−j∠A_(LComL11)) Γ_(S) ^(OInc2) [S^(OInc2)] Γ_(L) ^(OInc2)

When analyzing the signal-flow path for measured signal _(M)Dat₁₃, it is assumed that signal source 24 (FIG. 1) produces signal _(S)Tx₃₃ having a sinusoidal waveform (i.e., a frequency-domain signal) of frequency f_(S) and complex amplitude A_(S). For convenience, we include the effects of directional couplers 364, 366, 312, and 356 (FIG. 2), variable attenuator 302 (FIG. 2), and switches 304 (FIG. 2) and 106 (FIG. 3) in S-matrix [S^(Tx1)], that is modeled by S-parameters S₁₁ ^(Tx1), S₁₂ ^(Tx1), S₂₁ ^(Tx1), and S₂₂ ^(Tx1). Recall that [S^(Tx1)] also models the time-varying characteristics of Tx1 signal path 112 through high-frequency propagation medium 102. The error terms Γ_(S) ^(Tx1) and Γ_(L) ^(Tx1) model any mismatches at signal source 24 and at _(R)Dat₁₁ directional coupler 124 in first probe 100.

The discussion now assumes that all of components in first and second probes 100 and 200 (FIGS. 2 and 3) are ideal, i.e., they exhibit infinite return losses and the theoretical insertion losses (the time-varying effects of switch 106 has already been included in [S^(Tx1)]). An S-parameter {tilde over (S)}_(11M) represents the measured DUT 26 reflection coefficient, and it is assumed that the insertion loss along the path of signal Dat_(1J) through probe 100 (i.e., through directional couplers 122 and 124), the input pad, mixer, and output pad of mixing circuit 128, and through combiner 146) is characterized by a parameter P_(ODat11).

When signal _(R)Dat₁₁ enters data offset mixing circuit 128, it produces sum and difference frequencies, i.e., Equations (3) and (4), respectively. Since the sum frequency resides in the stopband of narrow-band crystal filter 454, only the difference-frequency signal will be included in the following discussion. In addition to the insertion loss of mixing circuit 128, the difference frequency signal _(O)Dat₁₁ will also experience a phase shift exp(−j∠A_(LDat11)), where ∠A_(LDat11) is the phase of local-oscillator signal _(L)Dat₁₁. This phase shift must be represented explicitly in FIG. 8 since it may vary over time, while the error term P_(ODat11) only accounts for static errors.

In order to simplify receiver Dat1 channel 404 analysis, the effects (return and insertion losses) of Dat1 channel 404 components (i.e., variable-gain amplifier 414, filter 424, IF mixer 434, pad 444, crystal filter 454, baseband mixer 464, amplifier 474, and lowpass filter 484) are included in the time-varying S-matrix [S^(Dat1)]. Here the signal undergoes two stages of down-conversion. However, this time we can include the phase shifts associated with IF signal source 438 and baseband signal source 468 in the time-varying S-matrix [S^(Dat1)] since these effects will be removed via signal-normalization techniques. The error terms Γ_(S) ^(ODat1) and Γ_(L) ^(ODat1) account for the mismatches at the output of probe 100 (FIG. 3) and the input to ADC 494 (FIG. 5), respectively.

Due to the fact that the insertion loss through probe 100 (FIG. 3) is large (i.e., −20 log₁₀|P_(ODat11)|>30 dB for all signal paths in probes 100 and 200), the signal that is reflected back into probe 100 by a mismatch in receiver Dat1 channel 404 will have negligible amplitude when it reaches transmitter 300 (FIG. 2). Thus, as shown in FIG. 8, the signal-flow path associated with the return signal in probe 100 is ignored, i.e., set to zero.

FIG. 9 shows a general signal-flow graph in accordance with a preferred embodiment of the present invention. The following discussion refers to FIGS. 2, 3, 5, 8, and 9.

The general signal-flow graph of Equation 11 may be used to investigate the effects of the transmitter and receiver circuits separately. The transmission coefficient T^(PathP) for this signal-flow graph can be written as

$\begin{matrix} \begin{matrix} {T^{PathP} = \frac{b_{2}^{PathP}}{b_{S}^{PathP}}} \\ {= {\frac{{\overset{\sim}{S}}_{21}^{PathP}}{\begin{matrix} {1 - {{\overset{\sim}{S}}_{11}^{PathP}\Gamma_{S}^{PathP}} - {{\overset{\sim}{S}}_{22}^{PathP}\Gamma_{L}^{PathP}} +} \\ {\left( {{{\overset{\sim}{S}}_{11}^{PathP}{\overset{\sim}{S}}_{22}^{PathP}} - {{\overset{\sim}{S}}_{21}^{PathP}{\overset{\sim}{S}}_{12}^{PathP}}} \right)\Gamma_{S}^{PathP}\Gamma_{L}^{PathP}} \end{matrix}}.}} \end{matrix} & (11) \end{matrix}$

Using this result, we find that the measured frequency-domain data signal _(M)Dat₁₃, which is part of composite signal _(B)Dat₁₃ sampled by ADC 494, can be written as _(M) Dat ₁₁ =A _(S) T ^(Tx1) {tilde over (S)} _(11M) P _(ODat11) exp(−j∠A _(LDat11))T ^(Dat1) G _(Dat1),  (12) where G_(Dat1) denotes the time-varying gain associated with the ADC 494, and T^(Tx1) and T_(Dat1) represent the time-varying transmission coefficients (Equation (11)) for the transmitter and data-receiver circuitry, respectively.

The result in Equation (12) is listed in the first row of Table 3. As discussed earlier, the signal-flow analyses for the other paths can be carried out in a similar manner, thus yielding the other results in Table 3 for the case of forward DUT signal measurements. The signal _(M)Inc₂₃ is not used for forward DUT signal measurements, but may be employed to monitor the time-varying load match, thereby avoiding the need for switch 206 and 50Ω loads 208 and 210 (FIG. 4).

TABLE 3 Signal Frequency-Domain Representation f_(PathP) _(M)Dat₁₃ A_(S)T^(Tx1){tilde over (S)}_(11M)P_(ODat11)exp(−j∠A_(LDat11))T^(Dat1)G_(Dat1) 510 Hz _(M)ComD₁₃ A_(S)T^(RefV1)P_(OComD11)exp(−j∠A_(LComD11))T^(Dat1)G_(Dat1) 810 Hz _(m)Inc₁₃ A_(S)T^(Tx1)P_(OInc11)exp(−j∠A_(LInc11))T^(Inc1)G_(Inc1) 700 Hz _(M)ComV₁₃ A_(S)T^(RefV1)P_(OComV11)exp(−j∠A_(LComV11))T^(Inc1)G_(Inc1) 860 Hz _(M)ComF₁₃ A_(S)T^(RefF1)P_(OComF11)exp(−j∠A_(LComF11))T^(Inc1)G_(Inc1) 310 Hz _(M)Link₂₃ AST^(RefF2)P_(Link22)exp(−j∠A_(LComL22)){tilde over (T)}^(Link)T^(Inc1)G_(Incl) 950 Hz _(M)Dat₂₃ A_(S)T^(Tx1){tilde over (S)}_(21M)P_(ODat22)exp(−j∠A_(LDat22))T^(Dat2)G_(Dat2) 680 Hz _(M)ComD₂₃ A_(S)T^(RefV2)P_(OComD22)exp(−j∠A_(LComD22))T^(Dat2)G_(Dat2) 690 Hz _(M)Inc₂₃ NA 980 Hz _(M)ComV₂₃ A_(S)T^(RefV2)P_(OComV22)exp(−j∠A_(LComV22))T^(Inc2)G_(Inc2) 590 Hz _(M)ComF₂₃ A_(S)T^(RefF2)P_(OComF22)exp(−j∠A_(LComF22))T^(Inc2)G_(Inc2) 250 Hz _(M)Link₁₃ A_(S)T^(RefF2)P_(Link11)exp(−j∠A_(LComL11)){tilde over (T)}^(Link)T^(Inc2)G_(Inc2) 990 Hz

Reverse DUT measurements can be made by switching switches 304 and 306 (FIG. 2), 106 (FIG. 3), and 206 (FIG. 4) to the reverse positions (opposite the positions shown in FIGS. 2, 3, and 4). This allows transmission signal Tx₃₃ to be directed to second probe 200. Due to the symmetry exhibited between forward and reverse DUT signal measurements, it is straightforward to determine the signal representations for the reverse case. These results are given in Table 4. For reverse signal measurements, signal _(M)Inc₁₁ remains unused.

TABLE 4 Signal Frequency-Domain Representation f_(PathP) _(M)Dat₁₃ A_(S)T^(Tx2){tilde over (S)}_(12M)P_(ODat11)exp(−j∠A_(LDat11))T^(Dat1)G_(Dat1) 510 Hz _(M)ComD₁₃ A_(S)T^(RefV1)P_(OComD11)exp(−j∠A_(LComD11))T^(Dat1)G_(Dat1) 810 Hz _(M)Inc₁₃ NA 700 Hz _(M)ComV₁₃ A_(S)T^(RefV1)P_(OComV11)exp(−j∠A_(LComV11))T^(Inc1)G_(Inc1) 860 Hz _(M)ComF₁₃ A_(S)T^(RefF1)P_(OComF11)exp(−j∠A_(LComF11)T^(Inc1)G_(Inc1) 310 Hz _(M)Link₂₃ A_(S)T^(RefF2)P_(Link22)exp(−j∠A_(LComL22)){tilde over (T)}^(Link)T^(Inc1)G_(Inc1) 950 Hz _(M)Dat₂₃ A_(S)T^(Tx2){tilde over (S)}_(22M)P_(ODat22)exp(−j∠A_(LDat22))T^(Dat2)G_(Dat2) 680 Hz _(M)ComD₂₃ A_(S)T^(RefV2)P_(OComD22)exp(−j∠A_(LComD22))T^(Dat2)G_(Dat2) 690 Hz _(M)Inc₂₃ A_(S)T^(Tx2)P_(OInc22)exp(−j∠A_(LInc22))T^(Inc2)G_(Inc2) 980 Hz _(M)ComV₂₃ A_(S)T^(RefV2)P_(OComV22)exp(−j∠A_(LComV22))T^(Inc2)G_(Inc2) 590 Hz _(M)ComF₂₃ A_(S)T^(RefF2)P_(OComF22)exp(−j∠A_(LComF22)T^(Inc2)G_(Inc2) 250 Hz _(M)Link₁₃ A_(S)T^(RefF2)P_(Link11)exp(−j∠A_(LComL1),){tilde over (T)}^(Link)T^(Inc2)G_(Inc2) 990 Hz

Reference to Tables 3 and 4 shows that the phases associated with local-oscillator signals _(L)Dat₁₁, _(L)ComD₁₁, _(L)Inc₁₁, _(L)ComV₁₁, _(L)ComF₁₁, and _(L)ComL₁₁ (FIG. 3), and _(L)Dat₂₂, _(L)ComD₂₂, _(L)Inc₂₂, _(L)ComV₂₂, _(L)ComF₂₂, and _(L)ComL₂₂ (FIG. 4) appear in composite receiving signals Inc_(1J), Dat_(1J), Dat_(2J), Inc₂₂, and Link_(PJ) (FIG. 5). Since these phase terms are not common to two or more signals, direct signal normalization techniques cannot be used to remove the associated time-varying errors. In order to address this problem, TES network analyzer 22 uses a separate ADC converter 492 in LOSum channel 102 to measure the relative phases between signals. These phases are then subtracted from the measured data by processor 500 during the signal normalization procedure. The frequency-domain signal representations associated with LOSum channel 102 are given in Table 5.

TABLE 5 Signal Frequency-Domain Representation f_(LPathP) _(M)Dat₁₃ |A_(LDat11)|exp(−j∠A_(LDat11))P_(LDat11)G^(LO) 99,800 Hz _(M)ComD₁₃ |A_(LComD11)|exp(−j∠A_(LComD11))P_(LComD11)G^(LO) 99,500 Hz _(M)Inc₁₃ |A_(LInc11)|exp(−j∠A_(LInc11))P_(LInc11)G^(LO) 99,610 Hz _(M)ComV₁₃ |A_(LComV11)|exp(−j∠A_(LComV11))P_(LComV11)G^(LO) 99,450 Hz _(M)ComF₁₃ |A_(LComF11)|exp(−j∠A_(LComF11))P_(LComF11)G^(LO) 100,000 Hz _(M)Link₂₃ |A_(LComL11)|exp(−j∠A_(LComL11))P_(LComL11)G^(LO) 99,360 Hz _(M)Dat₂₃ |A_(LDat22)|exp(−j∠A_(LDat22))P_(LDat22)G^(LO) 99,630 Hz _(M)ComD₂₃ |A_(LComD22)|exp(−j∠A_(LComD22))P_(LComD22)G^(LO) 99,620 Hz _(M)Inc₂₃ |A_(LInc22)|exp(−j∠A_(LInc22))P_(LInc22)G^(LO) 99,330 Hz _(M)ComV₂₃ |A_(LComV22)|exp(−j∠A_(LComV22))P_(LComV22)G^(LO) 99,720 Hz _(M)ComF₂₃ |A_(LComF22)|exp(−j∠A_(LComF22))P_(LComF22)G^(LO) 100,060 Hz _(M)Link₁₃ |A_(LComL22)|exp(−j∠A_(LComL22))P_(LComL22)G^(LO) 99,320 Hz

Analytical representations of the various signals that are extracted from the appropriate bins in the FFT after transforming the digital signals in TES network analyzer 22 having been obtained, signal normalization techniques may be used to remove both static and dynamic errors associated with the components through which both the data and normalization signals pass. Starting with forward DUT signal measurements, S-parameter S₁₁ involves the ratio of measured reflected DUT signal _(M)Dat₁₃ and measured incident DUT signal _(M)Inc₁₃. Reference to Table 3 shows that the ratio _(M)Dat₁₃:_(M)Inc₁₃ yields a result that is proportional to the desired S-parameter, but the result still contains dynamic errors associated with the phases of local-oscillator signals _(L)Dat_(1J) and _(L)Inc_(1J), the Dat1 and Inc1 receiver channels 404 and 403, and ADCs 494 and 493, respectively. However, this normalization process does remove any dynamic errors associated with signal source 24 (i.e., A_(S)), the port-1 transmitter circuitry (i.e., T^(Tx1)), and phase errors associated with IF signal source 438 and baseband signal source 468.

Further signal normalization is required to remove the remaining dynamic errors. If measured common signals _(M)ComD₁₃ and _(M)ComV₁₃ are included in the normalization, and their phases (via _(L)Sum₃₃ signal) are subtracted, then all of the dynamic error terms may be removed, i.e.,

$\begin{matrix} \begin{matrix} {S_{11M} = {\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{\exp\left\lbrack {j\left( {{\angle\; A_{LDat11}} + {\angle\; A_{LComV11}}} \right)} \right\rbrack}{\exp\left\lbrack {j\left( {{\angle\; A_{LComD11}} + {\angle\; A_{LInc11}}} \right)} \right\rbrack}}} \\ {= {{\overset{\sim}{S}}_{11M}\frac{P_{{MDat}\mspace{11mu} 13}}{P_{{MComD}\mspace{11mu} 13}}{\frac{P_{{MComV}\mspace{14mu} 13}}{P_{{MInc}\mspace{14mu} 13}}.}}} \end{matrix} & (13) \end{matrix}$

The static-error terms associated with the paths (i.e., P_(PathP)) in first probe 100 (FIG. 3) have been absorbed into the definition of S_(11M). As is shown hereinafter, these static errors are removed via static-error suppression techniques, thereby yielding the true S-parameter, S_(11A).

A similar procedure may be employed to find the second forward S-parameter, S₂₁. However, since this S-parameter involves the ratio of the transmitted wave _(M)Dat₂₃ and the incident wave _(M)Inc₁₃, the normalization signals on the two sides of DUT 26 must be linked together. This is accomplished via reciprocal link leg Link_(PJ) (i.e., Link channel 401 and Link signal path 140/240) and fixed-common signals ComF_(PP) in TES network analyzer 22. The overall transmission coefficient of reciprocal link leg Link_(PJ) is T^(Link). If the dynamic transmission coefficients [S^(Link1)] (FIG. 3) and [S^(Link2)] (FIG. 4) associated with the reciprocal link leg can be measured, then the result may be used to provide the necessary amplitude and phase references between first and second ports 30 and 32 of DUT 26. It may be shown that transmission coefficient T^(Link) may be accurately measured by taking the ratio of four signals as shown below:

$\begin{matrix} \begin{matrix} {T^{Link} = \sqrt{\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{\exp\left\lbrack {j\left( {{\angle\; A_{LComL22}} + {\angle\; A_{LComL11}}} \right)} \right\rbrack}{\exp\left\lbrack {j\left( {{\angle\; A_{LComF11}} + {\angle\; A_{LComF22}}} \right)} \right\rbrack}}} \\ {= {{\overset{\sim}{T}}^{Link}{\sqrt{\frac{P_{{MLink}\mspace{14mu} 23}}{P_{{MComF}\mspace{14mu} 13}}\frac{P_{{MLink}\mspace{14mu} 13}}{P_{{MComF}\mspace{14mu} 23}}}.}}} \end{matrix} & (14) \end{matrix}$

Once again, the static errors associated with the paths in first and second probes 100 and 200 are lumped into T^(Link). Note that all of the time-varying errors associated with the individual signals in Equation (14) are removed by the above normalization procedure (see Table 3).

Now that the link has been established between first and second ports 30 and 32, the dynamic errors can be removed from the measured DUT transmission coefficient using

$\begin{matrix} \begin{matrix} {S_{21M} = {\frac{1}{T^{Link}}\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}}} \\ {\frac{\exp\left\lbrack {j\left( {{\angle\; A_{LComL22}} + {\angle\; A_{LDat22}} + {\angle\; A_{LComV22}}} \right)} \right\rbrack}{\exp\left\lbrack {j\left( {{\angle\; A_{LInc11}} + {\angle\; A_{LComD22}} + {\angle\; A_{LComF22}}} \right)} \right\rbrack}} \\ {= {{\overset{\sim}{S}}_{21M}\frac{P_{{MDat}\mspace{14mu} 23}}{P_{{MComD}\mspace{14mu} 23}}\frac{P_{{MComV}\mspace{14mu} 23}}{P_{{MInc}\mspace{14mu} 13}}{\sqrt{\frac{P_{{MLink}\mspace{14mu} 23}}{P_{{MLink}\mspace{14mu} 13}}\frac{P_{{MComV}\mspace{14mu} 23}}{P_{{MComF}\mspace{14mu} 23}}}.}}} \end{matrix} & (15) \end{matrix}$

Since all of the signals that are employed in Equations (13), (14), and (15) are measured simultaneously by ADCs 492–497, both S-parameters may be measured simultaneously.

A similar procedure may be used to measure the reverse S-parameters. This time the signal representations in Table 4 are employed. The normalization equations for the reverse S-parameters are given below:

$\begin{matrix} \begin{matrix} {S_{22M} = {\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{\exp\left\lbrack {j\left( {{\angle\; A_{LDat22}} + {\angle\; A_{LComV22}}} \right)} \right\rbrack}{\exp\left\lbrack {j\left( {{\angle\; A_{LComD22}} + {\angle\; A_{LInc22}}} \right)} \right\rbrack}}} \\ {{= {{\overset{\sim}{S}}_{22M}\frac{P_{{MDat}\mspace{11mu} 23}}{P_{{MComD}\mspace{11mu} 23}}\frac{P_{{MComV}\mspace{14mu} 23}}{P_{{MInc}\mspace{14mu} 23}}}};} \end{matrix} & (16) \\ \begin{matrix} {T^{Link} = \sqrt{\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{\exp\left\lbrack {j\left( {{\angle\; A_{LComL22}} + {\angle\; A_{LCom11}}} \right)} \right\rbrack}{\exp\left\lbrack {j\left( {{\angle\; A_{LComF11}} + {\angle\; A_{LComF22}}} \right)} \right\rbrack}}} \\ {{= {{\overset{\sim}{T}}^{Link}\sqrt{\frac{P_{{MLink}\mspace{14mu} 23}}{P_{{MComF}\mspace{14mu} 13}}\frac{P_{{MLink}\mspace{14mu} 13}}{P_{{MComF}\mspace{14mu} 23}}}}};{and}} \end{matrix} & (17) \\ \begin{matrix} {S_{12M} = {\frac{1}{T^{Link}}\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}}} \\ {\frac{\exp\left\lbrack {j\left( {{\angle\; A_{LComL11}} + {\angle\; A_{LDat11}} + {\angle\; A_{LComV11}}} \right)} \right\rbrack}{\exp\left\lbrack {j\left( {{\angle\; A_{Linc22}} + {\angle\; A_{LComD11}} + {\angle\; A_{LComF11}}} \right)} \right\rbrack}.} \end{matrix} & (18) \end{matrix}$

Since signal-normalization techniques suppress all the static and dynamic linear errors associated with the components through which both the data and normalization signals pass, the only errors that remain after signal normalization are those in probes 100 and 200. In the analysis hereinbefore, these components were assumed to be ideal.

FIG. 10 shows a two-port 12-term error-correction model in accordance with a preferred embodiment of the present invention. The error-correction model of FIG. 10 accounts for directivity errors E_(DF) and E_(DR), reflection tracking errors E_(RF) and E_(RR), transmission tracking errors E_(TF) and E_(TR), source mismatch errors E_(SF) and E_(SR), load mismatch errors E_(LF) and E_(LR), and isolation errors E_(XF) and E_(XR).

Since the errors modeled by these systematic error terms must be repeatable, it is crucial to minimize the drift errors associated with all of the components in probes 100 and 200 (FIGS. 3 and 4). For this reason, the mixing circuits in probes 100 and 200 use double-balanced, ring-diode mixers, known to be very stable devices exhibiting minimal amplitude and phase drift. Furthermore, all other components in probes 100 and 200 are passive devices, and therefore do not contribute to significant signal amplitude and phase drift.

In applications that demand either high-accuracy or high-frequency measurements, where even a small amount of thermal expansion may lead to noticeable phase shift, probes 100 and 200 may be temperature stabilized. Temperature stabilization of probes 100 and 200 would be relatively easy since none of the components therein dissipate large amounts of heat and are all located in a relatively small volume. If temperature-stabilization is employed for probes 100 and 200, then it is not necessary to closely regulate ambient room temperature. Either way, it may be assumed that the circuit components in probes 100 and 200 produce only repeatable systematic errors, and can therefore be modeled by the previously discussed systematic-error terms.

Application of signal-flow analysis techniques to the error-correction model in FIG. 10 yields the expressions for the “actual” S-parameters for DUT 26 listed below as

$\begin{matrix} {{S_{11A} = \frac{\begin{matrix} {{\left( \frac{S_{11M} - E_{DF}}{E_{RF}} \right)\left\lbrack {1 + {\left( \frac{S_{22M} - E_{DR}}{E_{RR}} \right)E_{SR}}} \right\rbrack} -} \\ {\left( \frac{S_{21M} - E_{XF}}{E_{TF}} \right)\left( \frac{S_{12M} - E_{XR}}{E_{TR}} \right)E_{LF}} \end{matrix}}{\begin{matrix} {{\left\lbrack {1 + {\left( \frac{S_{11M} - E_{DF}}{E_{RF}} \right)E_{SF}}} \right\rbrack\left\lbrack {1 + {\left( \frac{S_{22M} - E_{DR}}{E_{RR}} \right)E_{SR}}} \right\rbrack} -} \\ {\left( \frac{S_{21M} - E_{XF}}{E_{TF}} \right)\left( \frac{S_{12M} - E_{XR}}{E_{TR}} \right)E_{LF}E_{LR}} \end{matrix}}},} & (19) \\ {{S_{21A} = \frac{\left\lbrack {1 + {\left( \frac{S_{22M} - E_{DR}}{E_{RR}} \right)\left( {E_{SR} - E_{LF}} \right)}} \right\rbrack\left( \frac{S_{21M} - E_{XF}}{E_{TF}} \right)}{\begin{matrix} {{\left\lbrack {1 + {\left( \frac{S_{11M} - E_{DF}}{E_{RF}} \right)E_{SF}}} \right\rbrack\left\lbrack {1 + {\left( \frac{S_{22M} - E_{DR}}{E_{RR}} \right)E_{SR}}} \right\rbrack} -} \\ {\left( \frac{S_{21M} - E_{XF}}{E_{TF}} \right)\left( \frac{S_{12M} - E_{XR}}{E_{TR}} \right)E_{LF}E_{LR}} \end{matrix}}},} & (20) \\ \begin{matrix} {{S_{12A} = \frac{\left\lbrack {1 + {\left( \frac{S_{11M} - E_{DF}}{E_{RF}} \right)\left( {E_{SF} - E_{LR}} \right)}} \right\rbrack\left( \frac{S_{12M} - E_{XR}}{E_{TR}} \right)}{\begin{matrix} {{\left\lbrack {1 + {\left( \frac{S_{11M} - E_{DF}}{E_{RF}} \right)E_{SF}}} \right\rbrack\left\lbrack {1 + {\left( \frac{S_{22M} - E_{DR}}{E_{RR}} \right)E_{SR}}} \right\rbrack} -} \\ {\left( \frac{S_{21M} - E_{XF}}{E_{TF}} \right)\left( \frac{S_{12M} - E_{XR}}{E_{TR}} \right)E_{LF}E_{LR}} \end{matrix}}},} \\ {and} \end{matrix} & (21) \\ {S_{22A} = {\frac{\begin{matrix} {{\left( \frac{S_{22M} - E_{DR}}{E_{RR}} \right)\left\lbrack {1 + {\left( \frac{S_{11M} - E_{DF}}{E_{RF}} \right)E_{SF}}} \right\rbrack} -} \\ {\left( \frac{S_{21M} - E_{XF}}{E_{TF}} \right)\left( \frac{S_{12M} - E_{XR}}{E_{TR}} \right)E_{LR}} \end{matrix}}{\begin{matrix} {{\left\lbrack {1 + {\left( \frac{S_{11M} - E_{DF}}{E_{RF}} \right)E_{SF}}} \right\rbrack\left\lbrack {1 + {\left( \frac{S_{22M} - E_{DR}}{E_{RR}} \right)E_{SR}}} \right\rbrack} -} \\ {\left( \frac{S_{21M} - E_{XF}}{E_{TF}} \right)\left( \frac{S_{12M} - E_{XR}}{E_{TR}} \right)E_{LF}E_{LR}} \end{matrix}}.}} & (22) \end{matrix}$

If a SOLT calibration technique is employed, then the static error coefficients can be computed using

$\begin{matrix} {{E_{DF} = S_{11M}^{L}},} & (23) \\ {{E_{DR} = S_{22M}^{L}},} & (24) \\ {{E_{RF} = \frac{{2\left\lbrack {S_{11M}^{O} - E_{DF}} \right\rbrack}\left\lbrack {S_{11M}^{S} - E_{DF}} \right\rbrack}{S_{11M}^{S} - S_{11M}^{O}}},} & (25) \\ {{E_{RR} = \frac{{2\left\lbrack {S_{22M}^{O} - E_{DR}} \right\rbrack}\left\lbrack {S_{22M}^{S} - E_{DR}} \right\rbrack}{S_{22M}^{S} - S_{22M}^{O}}},} & (26) \\ {{E_{SF} = \frac{{2E_{DF}} - S_{11M}^{O} - S_{11M}^{S}}{S_{11M}^{S} - S_{11M}^{O}}},} & (27) \\ {{E_{SR} = \frac{{2E_{DR}} - S_{22M}^{O} - S_{22M}^{S}}{S_{22M}^{S} - S_{22M}^{O}}},} & (28) \\ {{E_{XF} = S_{21M}^{I}},} & (29) \\ {{E_{XR} = S_{12M}^{I}},} & (30) \\ {{E_{TF} = \frac{E_{RF}\left( {S_{21M}^{T} - E_{XF}} \right)}{{E_{SF}\left( {S_{11M}^{T} - E_{DF}} \right)} + E_{RF}}},} & (31) \\ {{E_{TR} = \frac{E_{RR}\left( {S_{12M}^{T} - E_{XR}} \right)}{{E_{SR}\left( {S_{22M}^{T} - E_{DR}} \right)} + E_{RR}}},} & (32) \\ {{E_{LF} = {\frac{1}{E_{SF}} - \frac{E_{RF}}{\left\lbrack {{E_{SF}\left( {S_{11M}^{T} - E_{DF}} \right)} + E_{RF}} \right\rbrack E_{SF}}}},{and}} & (33) \\ {E_{LR} = {\frac{1}{E_{SR}} - {\frac{E_{RR}}{\left\lbrack {{E_{SR}\left( {S_{22M}^{T} - E_{DR}} \right)} + E_{RR}} \right\rbrack E_{SR}}.}}} & (34) \end{matrix}$

However, before Equations (19)–(22) may be employed, the error coefficients must be calculated by making a set of measurements on a set of known loads (i.e., precise standards). A sufficient number of precise standards must be used in order to back-out the various error coefficients in the error model. The most common technique involves measurement of a Short circuit, an Open circuit, a matched Load, and a Through (SOLT). If we denote the measured S-parameters for a given DUT by S_(piM) ^(d), where i=1,2 denotes the transmitter port, p=1,2 denotes the receiver port, and d=S,O,L,T,I denotes the type of DUT employed in the calibration, then the expressions Equations (23)–(34) for the systematic-error terms can be calculated from the measured S-parameters. The case d=I indicates that the proper DUTs have been connected for the isolation measurements.

Note that the results that are computed using Equations (13), (15), (16), and (18) are substituted for the measured S-parameters in Equations (19)–(34).

Those skilled in the art will appreciate that while a full two-port error model has been employed in the TES-technique proof-of-principle exemplary embodiment discussed herein, simpler error models (e.g., one-port reflection or 2-port transmission response error models) may also have been employed in the TES technique. The use of other error models does not depart from the spirit of the present invention.

Random errors may consist of thermal noise, shot noise, and electromagnetic interference. The TES incorporates random-error suppression techniques such as stacking or signal averaging, filtering, oversampling, lock-in amplifiers, and synchronous detection. In all of these approaches, substantial random-noise reduction requires long integration times. With the TES technique, system response drift is eliminated. This allows the use of virtually unlimited integration times and therefore unlimited reduction in random noise errors.

Phase noise is a performance-limiting factor in applications where a weak signal must be detected in the presence of a stronger, interfering signal. In TES network analyzer 22, multiple signals reside at different frequencies and varying signal levels, and are simultaneously detected on a single receiver channel. TES network analyzer 22 minimizes the effects of phase noise interference by maintaining the levels of signals _(M)Inc₁₃, _(M)ComV₁₃, _(M)ComF₁₃, and _(M)ComL₂₃ (_(M)Link₂₃), signals _(M)Dat₁₃ and _(M)ComD₁₃, signals _(M)Dat₂₃ and _(M)ComD₁₃, and signals _(M)Inc₂₃, _(M)ComV₂₃, _(M)ComF₂₃, and _(M)ComL₁₃ (_(M)Link₁₃) within composite signals Inc_(1J), Dat_(1J), Dat_(2J), and Inc_(2J), respectively, within ±15 dB of each other. However, phase noise can still present problems in high-frequency applications where the source has a phase-noise level that is greater than −75 dBc at the location of the next frequency bin having a valid signal (i.e., after Fourier transformation). In such cases, a phase-noise error suppression by subtraction technique may be employed.

The basic idea behind the phase-noise error suppression by subtraction technique is to measure the complex amplitude of the phase noise relative to source signal _(S)Tx₃₃ in dedicated PN channel 407 containing a single signal. This information is then used to extract the phase-noise contamination in other channels 403–406 containing multiple signals. It is assumed that the differences between the time delays associated with signal propagation along channels 403–407 are small in comparison to the data acquisition time of ADCs 493–497. Furthermore, it is assumed that all signals are simultaneously measured by ADCs 493–497. With these assumptions, the phase-noise suppression by subtraction technique works since all signals are derived from a common high-frequency signal Tx₃₃ whose source 24 is the cause of the phase-noise contamination. Down-conversion sources 438 and 468 may also contribute appreciable amounts of phase noise to the system. However, phase noise contributed by down-conversion sources 438 and 468 is common to all signals and may be monitored and removed. Signals from local-oscillator module 354, which are not common to all measured signals, contribute negligible levels of phase-noise contamination, and may therefore be ignored.

Referring to FIGS. 2 and 5, the effects of down-conversion and filtering in PN channel 407 may be readily understood. Extracted phase-noise signal _(X)PN₃₃ is extracted from source signal _(S)Tx₃₃ by coupler 356. Signal _(X)PN₃₃ resides at substantially predetermined frequency f_(S) of source signal _(S)Tx₃₃. Mixing circuit 358 mixes signal _(X)PN₃₃ with local-oscillator phase-noise offset signal _(L)PN₃₃ having a frequency f_(LO). This produces the upper and lower sidebands at frequencies f_(S)+f_(LO) and f_(S)−f_(LO) (see Equations (3) and (4)).

It is the lower sideband having a frequency of f_(S)−f_(LO) that is desired (the upper sideband will be filtered out). It is assumed signal source 24 produces −80 dBc/Hz phase noise 10 kHz from predetermined frequency f_(S), and −110 dBc/Hz phase noise 100 kHz from f_(S). Thus, using offset local-oscillator signals with frequencies of approximately 100 kHz assures that the phase noise contamination from the upper sideband will be at least −110 dBc/Hz at the location of the desired lower sideband.

IF down-conversion mixer 437 mixes the upper and lower sideband signals with IF source signal _(S)IF₃₃ having a frequency of f_(IF). This produces four distinct signals having frequencies of f_(S)−f_(IF)−f_(LO), f_(S)−f_(IF)+f_(LO), f_(S)+f_(IF)−f_(LO), and f_(S)+f_(IF)+f_(LO) (see Equations (5)–(8)). High-Q crystal bandpass filter 457 has a 2 kHz bandwidth and passes the desired signal of frequency f_(S)−f_(IF)−f_(LO), while rejecting the other unwanted signals including the phase noise that wraps around DC.

The relative level of the phase noise that is measured in PN channel 407 at a frequency offset f_(Δ) from the carrier frequency is defined as

$\begin{matrix} {{\Phi_{M}\left( f_{\Delta} \right)} = {\frac{{Sig}\;{\Phi_{M}\left( {f_{carrier} - f_{\Delta}} \right)}}{{Sig}\;{\Phi_{M}\left( f_{carrier} \right)}}.}} & (35) \end{matrix}$

In the above expression, SigΦ_(M)(f) denotes the signal that is measured on PN channel 407 at the frequency f. Note that the normalization in Equation (35) removes the amplitude and phase dependencies associated with signal source 24, IF signal source 438, baseband signal source 468, and local-oscillator module 354. One of the other receiver channels 403–406 is chosen that contains multiple simultaneous frequency-offset signals, i.e., K signals. The frequency-domain representation for the J-th frequency-offset signal, measured at its carrier frequency f_(J), can be represented analytically by

$\begin{matrix} {{{{Sig}\;{J_{M}\left( f_{J} \right)}} = {{{Sig}\;{J_{A}\left( f_{J} \right)}} + {\sum\limits_{\underset{K \neq J}{K = 1}}^{K}{C_{J}^{K}{\Phi_{M}\left( {f_{K} - f_{J}} \right)}{Sig}\;{K_{A}\left( f_{K} \right)}}}}};} & (36) \end{matrix}$ where J=1,2, . . . ,K, and where SigJ_(A)(f_(J)) denotes the actual (i.e., uncontaminated) value of the J-th signal at the frequency f_(J), SigK_(A)(f_(K)) denotes the actual value of the K-th (for K≠J) signal at its carrier frequency f_(K), and the constants C_(J) ^(K) model the static errors that result from the nonideal and nonidentical components that are used in this receiver channel and PN channel 407. The above equation can be written in matrix form as

$\begin{matrix} {\begin{bmatrix} {{Sig}\; 1_{M}\left( f_{1} \right)} \\ {{Sig}\; 2_{M}\left( f_{2} \right)} \\ \vdots \\ {{Sig}\;{K_{M}\left( f_{K} \right)}} \end{bmatrix} = {\quad{\begin{bmatrix} 1 & {C_{1}^{2}{\Phi_{M}\left( {f_{2\;} - f_{1}} \right)}} & \ldots & {C_{1}^{K}{\Phi_{M}\left( {f_{K} - f_{1}} \right)}} \\ {C_{2}^{1}{\Phi_{M}\left( {f_{1} - f_{2}} \right)}} & 1 & \ldots & {C_{2}^{K}{\Phi_{M}\left( {f_{K} - f_{2}} \right)}} \\ \vdots & \vdots & ⋰ & \vdots \\ {C_{K}^{1}{\Phi_{M}\left( {f_{1} - f_{K}} \right)}} & {C_{K}^{2}{\Phi_{M}\left( {f_{2} - f_{K}} \right)}} & \ldots & 1 \end{bmatrix}{\quad{\begin{bmatrix} {{Sig}\; 1_{A}\left( f_{1} \right)} \\ {{Sig}\; 2_{A}\left( f_{2} \right)} \\ \vdots \\ {{Sig}\;{K_{A}\left( f_{K} \right)}} \end{bmatrix}.}}}}} & (37) \end{matrix}$

The desired actual signal levels can be determined from the measured phase-noise contaminated signals by solving this system of simultaneous equations, i.e.,

$\begin{matrix} {\begin{bmatrix} {{Sig1}_{A}\left( f_{1} \right)} \\ {{Sig2}_{A}\left( f_{2} \right)} \\ \vdots \\ {{SigK}_{A}\left( f_{K} \right)} \end{bmatrix} = {{\begin{bmatrix} 1 & {C_{1}^{2}{\Phi_{M}\left( {f_{2} - f_{1}} \right)}} & \cdots & {C_{1}^{K}{\Phi_{M}\left( {f_{K} - f_{1}} \right)}} \\ {C_{2}^{1}{\Phi_{M}\left( {f_{1} - f_{2}} \right)}} & 1 & \cdots & {C_{2}^{K}{\Phi_{M}\left( {f_{K} - f_{2}} \right)}} \\ \vdots & \vdots & ⋰ & \vdots \\ {C_{K}^{1}{\Phi_{M}\left( {f_{1} - f_{K}} \right)}} & {C_{K}^{2}{\Phi_{M}\left( {f_{2} - f_{K}} \right)}} & \cdots & 1 \end{bmatrix}^{- 1}\left\lbrack \begin{matrix} {{Sig1}_{M}\left( f_{1} \right)} \\ {{Sig2}_{M}\left( f_{2} \right)} \\ \vdots \\ {{SigK}_{M}\left( f_{K} \right)} \end{matrix} \right\rbrack}.}} & (38) \end{matrix}$

Before this matrix equation can be employed, the static calibration constants, C_(J) ^(K), must be found. These constants are obtained during a static calibration procedure that is performed prior to taking measurements. This phase-noise calibration procedure is carried out before the data have been collected for the 12-term static calibration procedure. In order to extract the calibration constants, each of the signals are collected separately in the absence of the other interfering signals. If all of the signals, with the exception of the U-th signal, are turned off in the receiver channel, then Equation (36) can be used to show that the ratio of the phase noise at the location of the V-th signal and the carrier signal can be written as

$\begin{matrix} \begin{matrix} {{Rat}_{V/U}^{Uonly} = \left\lbrack \frac{{SigV}_{M}\left( f_{V} \right)}{{SigU}_{M}\left( f_{U} \right)} \right\rbrack^{Uonly}} \\ {= {C_{V}^{U}{{\Phi_{M}\left( {f_{U} - f_{V}} \right)}.}}} \end{matrix} & (39) \end{matrix}$

The desired calibration coefficients can then be obtained via

$\begin{matrix} {{C_{V}^{U} = \frac{{Rat}_{V/U}^{Uonly}}{\Phi_{M}\left( {f_{U} - f_{V}} \right)}},} & (40) \end{matrix}$ where Φ_(M)(f_(U)−f_(V)) is measured in PN channel 407 (see Equation (35)). Note that signals Dat_(PJ) and Inc_(PJ) may be set to zero by opening switches 106 (FIG. 3) and 206 (FIG. 4). Signals ComD_(PP) and ComV_(PP) (FIGS. 3 and 4) may be removed using variable attenuators 322 and 324 (FIG. 2). In a similar manner, signals ComF_(PP) and Link_(PJ) (FIGS. 3 and 4) may be removed using variable attenuators 318 and 320 (FIG. 2). Unfortunately, when one of these attenuators is set to a large value, both signals ComF_(PP) and Link_(PJ) are removed. In order to only turn off one of these signals during the calibration process, the appropriate local-oscillator source may be turned off, thereby turning the corresponding mixing circuit into an attenuator. As long as the calibration procedure is carried out relatively quickly, this lack of applied local-oscillator signal will not result in measurable temperature or signal drift in the mixing circuit.

In order to accurately measure the calibration coefficients via Equation (40), it is necessary to accurately measure both

Rat_(V/U)^(Uonly) and Φ_(M)(f_(U)−f_(V)). Since the phase-noise sidebands are typically small relative to the size of the carrier, it is difficult to measure these signals accurately. In order to improve the accuracy of these measurements, mixer 362 can be used to simulate higher level phase noise side-bands during the calibration procedure. The local-oscillator signal frequency for mixer 362 is set to f_(LΦCal)=f_(U)−θ_(V) during the calibration procedure.

The TES technique involves a multilevel calibration technique that consists of static, dynamic, phase-noise, and random error suppression techniques. After the time-domain signals are measured by ADCs 492–497, processor 500 (FIG. 7) uses FFT routine 510 to find the complex amplitudes for frequency-domain signals _(M)Dat₁₃, _(M)ComD₁₃, _(M)Inc₁₃, _(M)ComV₁₃, _(M)ComF₁₃, _(M)ComL₂₃ (_(M)Link₂₃), _(M)Dat₂₃, _(M)ComD₂₃, _(M)Inc₂₃, _(M)ComV₂₃, _(M)ComF₂₃, and _(M)ComL₂₃ (_(M)Link₁₃), as listed in Table 3. In order to remove the phase-noise contamination from these signals, the phase-noise suppression by subtraction technique (see Equation (38)) is applied to signals _(M)Inc₁₃, _(M)ComV₁₃, _(M)ComF₁₃, and _(M)ComL₂₃ (_(M)Link₂₃), signals _(M)Dat₁₃ and _(M)ComD₁₃, signals _(M)Dat₂₃ and _(M)ComD₁₃, and signals _(M)Inc₂₃, _(M)ComV₂₃, _(M)ComF₂₃, and _(M)ComL₁₃ (_(M)Link₁₃) within composite signals Inc_(1J), Dat_(1J), Dat_(2J), and Inc_(2J), respectively, in receiver channels 403, 404, 405, and 406. The signal and its associated phase-noise sidebands appearing in PN channel 407 are used to compute the relative phase noise (see Equation (35)) that is required in Equation (38). Once computed using Equation (38), then the phase-noise-corrected, actual signals are inserted into the TES signal normalization equations (13)–(18), which remove the dynamic errors and some of the static errors, thereby yielding the measured S-parameters. The remaining static errors associated with the components in probes 100 and 200 are suppressed by substituting the calculated, measured S-parameters into equations (19)–(22), thus producing the actual S-parameters for DUT 26.

FIG. 11 shows a simplified schematic diagram of a reciprocal amplification circuit 600 for Link signal path 140/240 of TES network analyzer 22, and FIG. 12 shows a chart depicting relative frequencies for signals in reciprocal amplification circuit 600 in accordance with an alternative preferred embodiment of the present invention. The following discussion refers to FIGS. 11 and 12.

As is discussed hereinbefore, signal normalization techniques can be used to suppress static and dynamic errors associated with the components in first probe 100 (FIG. 3), second probe 200 (FIG. 4), transmitter 300 (FIG. 2), and receiver 400 (FIG. 5) through which both the data and normalization signals flow. In TES network analyzer 26, a reciprocal link leg (Link signal path 140/240 in FIGS. 3, 4, and 5) was employed to achieve accurate phase and magnitude references between ports 100 and 200. Two signals Link_(1J) and Link_(2J) (where J=1,2,3), which exist at different frequencies, simultaneously flow in opposite directions through the reciprocal link leg. The amplitude and phase references between probes 100 and 200 (i.e., the transmission coefficient T^(Link) on the link leg) are measured by combining link signals Link_(1J) and Link_(2J) with fixed-reference signals, as shown in Equation (14).

The pads for signals Link₁₁, Link₂₁, and Link_(P1), i.e. the output pad of mixing circuit 130, pad 162, and pad 160, respectively, serve to isolate first probe 100 from time-variations in the reflection properties associated with Link1 signal path 140 (FIG. 3). Similarly, the pads for signals Link₂₂, Link₁₂, and Link_(P2), i.e. the output pad of mixing circuit 230, pad 262, and pad 260, respectively, serve to isolate second probe 200 from time-variations in the reflection properties associated with Link2 signal path 240 (FIG. 4). Because of this, the Link_(PJ) (where P=1,2 and J=P,3) signal levels at the outputs of four-way combiners 142 and 242 is −62 dBm. When choosing the levels for the signals that are combined in the Inc1 and Inc2 receiver channels 403 and 406, this relatively low level for link signals Link_(PJ) places a limit on how large the other signal levels in channels 403 and 406 may be. This may cause a problem at high frequencies where link signals Link_(PJ) would suffer much greater cable losses. In this instance, it would be desirable to amplify link signals Link_(PJ).

The signal normalization procedure depends on having a reciprocal link leg, i.e., bidirectional Link signal path 140/240. Therefore, conventional amplifiers, being unidirectional, cannot be placed directly in Link signal path 140/240. Reciprocal-amplification circuit 600 of FIG. 11 may be employed. Note that this circuit will receive link signals Link_(1J) and Link_(2J) from probes 100 and 200, amplify them, and send them to probes 200 and 100, respectively, even though they are simultaneously propagating over Link signal path 140/240 in opposite directions. The presence of both signals Link_(1J) and Link_(2J) at both probes 100 and 200 will not pose difficulties since each signal Link_(1J) and Link_(2J) is used only in the probe 100 or 200 in which it was not generated, i.e., Link_(1J), being generated in first probe 100 (as signal Link₁₁), is only used in second probe 200 (as signal Link₁₂).

Referring to FIG. 11, reciprocal amplification circuit 600 is interposed within Link signal path 140/240 in receiver 400, dividing Link signal path 140/240 into a port-1 part (Link1 signal path 140) and a port-2 part (Link2 signal path 240).

To understand the operation of reciprocal-amplification circuit 600, it is necessary to follow the progress of signals Link_(1J) and Link_(2J) through the circuit. To facilitate this, “measurement” points 601, 602, 603, 604, 605, 606, 607, 608, 609, 610, 611, 612, 613, 614, 615, 616, 617, 618, 619, 620, 621, 622, 623, 624, 625, 626, 627, 628, 629, 630, 631, 633, 634, 635, and 636 have been placed in circuit 600, with points 606 and 631 serving as i/o ports for circuit 600. Also, an abbreviated nomenclature will be used to simplify this discussion.

Signals Link₁₁ and Link₂₂ have different frequencies. In the discussion of reciprocal-amplification circuit 600, these frequencies shall be referred to as f_(L1) and f_(L2), respectively, where f_(L1)≠f_(L2). Frequencies f_(L1) and f_(L2) are each derived from and source signal _(S)Tx₃₃ (FIGS. 2, 3, and 4) and are proximate source signal _(S)Tx₃₃ in frequency, i.e., frequencies f_(L1) and f_(L2) are nearly equal to predetermined frequency f_(S). Since, in the exemplary embodiment, source signal _(S)Tx₃₃ is a swept signal ranging in frequency from approximately 45 MHz to approximately 20 GHz, f_(L1) and f_(L2), too, range from approximately 45 MHz to 20 GHz.

Reciprocal-amplification circuit 600 uses filters to inhibit oscillation. Because of this, circuit 600 operates at fixed frequencies only. Therefore, frequencies f_(L1) and f_(L2) must be converted to known fixed frequencies for processing by circuit 600. This is accomplished via bidirectional hybrid mixers 640 and 642. For signals flowing towards circuit 600, mixers 640 and 642 function as down-conversion mixers. Conversely, for signals flowing away from circuit 600, mixer 640 and 642 function as up-conversion mixers.

Mixers 640 and 642 perform down-conversions and up-conversions by mixing frequencies f_(L1) and f_(L2) with a swept intermediate frequency source. In the exemplary embodiment, this swept frequency source is IF source signal _(S)IF₃₃, since, as described hereinbefore, signal _(S)IF₃₃ is a swept signal synchronized with source signal _(S)Tx₃₃, hence with link signals Link₁₁ and Link₂₂. As discussed in more detail below, the result of mixing frequencies f_(L1) and f_(L2) with IF source signal _(S)IF₃₃ is the frequencies f₁ and f₂, respectively.

Those skilled in the art will appreciate that the use of IF source signal _(S)IF₃₃ is not a requirement of the present invention, and that the use of other signals for down-conversion and up conversion does not depart from the spirit of the present invention.

Substantially midway between fixed frequencies f₁ and f₂ is a third frequency, f₀. The difference between f₁ and f₀ constitute a frequency interval. The scale of the chart in FIG. 12 is calibrated in such intervals. All frequencies in this discussion shall be presumed to be substantially equal to integral frequencies on this scale. That is, all frequencies shall be presumed to fall proximate the interval markers of FIG. 12.

The following discussion traces signal Link_(1J) through reciprocal-amplification circuit 600 in conjunction with Table 6. Table 6 depicts the frequencies of signal Link_(1J) at each of points 601, and 603–635 throughout reciprocal amplification circuit 600 (points 602 and 636 are not used for signal Link_(1J)).

TABLE 6 Points Frequencies of Link_(1J) 601 603 f_(L1) 604 f_(IF) 605 606 f_(L1), f_(IF) − f_(L1), f_(IF), f_(IF) + f_(L1) 607 608 609 610 611 f₁, i.e., f_(IF) − f_(L) ₁ 612 f_(A) 613 614 f_(A), f₁ − f_(A), f₁, f₁ + f_(A) 615 616 617 618 619 f₁ + f_(A) 620 621 622 623 624 f_(B) 625 626 f_(B), f₁ + f_(A) − f_(B), f₁ + f_(A), f₁ + f_(A) + f_(B) 627 628 629 630 631 f₁ + f_(A) − f_(B) 632 633 f_(IF) 634 635 f₁ + f_(A) − f_(B), f_(IF) − f₁ + f_(A) − f_(B), f_(IF), f_(IF) + f₁ + f_(A) − f_(B) Result f_(IF) + f₁ + f_(A) − f_(B)≈f_(L1)

Signal Link_(1J) originates in first probe 100 with frequency f_(L1), higher- and lower-frequency components.

Signal Link_(1J) propagates from point 601, through combiner/splitter 138 in first probe 100, through point 603, and into bidirectional hybrid mixer 640. Mixer 640 mixes signal Link_(1J) at frequency f_(L1) from point 603 with signal _(S)IF₃₃ at frequency f_(IF) from point 604 to produce signal Link_(1J) at frequencies f_(L1), f_(IF)−f_(L1), f_(IF), and f_(IF)+f_(L1) (plus higher-order components, which may be ignored) at point 605.

From point 605, signal Link_(1J) propagates over Link1 signal path 140 through high-frequency propagation medium 102 and into a low-high subcircuit 650 at point 606. Reciprocal-amplification circuit 600 is divided into subcircuits 650 and 652. Subcircuits 650 and 652 are identified as low-high subcircuit 650 (i.e., the subcircuit loop contains a low-center filter before a high-center filter) and high-low subcircuit 652 (i.e., the subcircuit loop contains a high-center filter before a low-center filter). Subcircuits 650 and 652 are coupled together by a propagation medium 654 having a LinkP signal path 656.

From point 606, signal Link_(1J) propagates through bandpass filter 658 to point 607. Bandpass filter 658 has a passband slightly greater than two intervals and centered on frequency f₀. That is, bandpass filter 658 is configured to pass frequency f₀, to pass frequency f₁ located one interval below frequency f₀, to pass frequency f₂ located one interval above frequency f₀, and to attenuate all other frequencies located two or more intervals either below or above frequency f₀. Since f₁=f_(1H)−f_(L1), bandpass filter 658 passes signal Link_(1J) at frequency f₁ while attenuating the other components.

From point 607, signal Link_(1J) propagates through combiner/splitter 660, through point 608, through combiner/splitter 662, through point 609, and into amplifier 664. Amplifier 664 serves not only to amplify signal Link_(1J), but also as a unidirectional isolator, i.e., to prohibit reverse propagation of signals from point 610 back to point 609.

From amplifier 664, signal Link_(1J) propagates through point 610 and into low-center bandpass filter 666. Bandpass filter 666 has a passband slightly greater than two intervals and centered on a frequency two intervals below frequency f₀, (i.e., low-center relative to f₀). That is, bandpass filter 650 is configured to pass frequencies located at three, two, and one intervals below frequency f₀, and to attenuate frequencies located four or more intervals below frequency f₀, or at or above frequency f₀. Since frequency f₁ is one interval below frequency f₀, bandpass filter 650 passes signal Link_(1J).

From bandpass filter 666, signal Link_(1J) propagates through point 611 and into mixer 668. Mixer 668 mixes signal Link_(1J) at frequency f₁ from point 611 with a signal (not shown) at frequency f_(A) from point 612 to produce signal Link_(1J) at frequencies f_(A), f₁−f_(A), f₁, and f₁+f_(A) (plus higher-order components, which may be ignored) at point 613. Frequency f_(A) is substantially equal to three intervals. Therefore, since frequency f₁ is substantially one interval below f₀, frequency f₁−f_(A) is substantially four intervals below f₀ and frequency f₁+f_(A) is substantially two intervals above f₀.

From point 613, signal Link_(1J) propagates into an amplifier 670. Like amplifier 664, amplifier 670 serves not only to amplify signal Link_(1J), but also as a unidirectional isolator, i.e., to prohibit reverse propagation of signals from point 614 back to point 613.

From amplifier 670, signal Link_(1J) propagates through point 614 and into high-center bandpass filter 672. Bandpass filter 672 has a passband slightly greater than two intervals and centered on a frequency two intervals above frequency f₀ (i.e., high-center relative to f₀). That is, bandpass filter 672 is configured to pass frequencies located at one, two, and three intervals above frequency f₀, and to attenuate frequencies located four or more intervals above frequency f₀, or at or below frequency f₀. Since frequency f₁+f_(A) is two intervals above frequency f₀ while frequencies f_(A), f₁−f_(A), and f₁ are all below frequency f₀, frequency f₁+f_(A) is passed and frequencies f_(A), f₁−f_(A), and f₁ are all attenuated.

From filter 672, signal Link_(1J) propagates through point 615 and into amplifier 674. Like amplifiers 664 and 670, amplifier 674 serves not only to amplify signal Link_(1J), but also as a unidirectional isolator, i.e., to prohibit reverse propagation of signals from point 616 back to point 615.

From amplifier 674, signal Link_(1J) propagates through point 616 and into combiner/splitter 662. The circuit from combiner/splitter 662, through amplifier 664, through filter 666, through mixer 668, through amplifier 670, through filter 672, through amplifier 674, and back to combiner/splitter 662 forms a loop. Filters 666 and 672, in conjunction with mixer 668, inhibit this loop from oscillating.

From combiner/splitter 662, signal Link_(1J) propagates through point 617 (coincident with point 608), through combiner splitter 660, through point 618, through propagation medium 654 via LinkP signal path 656, into high-low subcircuit 652, through point 619, through a combiner/splitter 676, through point 620, through a combiner/splitter 678, through point 621, and into an amplifier 680. From amplifier 680, signal Link_(1J) propagates through point 622 and into a bandpass filter 682. Like bandpass filter 672 discussed hereinbefore, bandpass filter 682 has a passband slightly greater than two intervals and centered on a frequency two intervals above frequency f₀. Since frequency f₁+f_(A) is two intervals above frequency f₀, frequency f₁+f_(A) is passed.

From bandpass filter 682, signal Link_(1J) propagates through point 623 and into a mixer 684. Within mixer 684, signal Link_(1J) at frequency f₁+f_(A) from point 623 is mixed with a signal (not shown) at frequency f_(B) from point 622 to produce signal Link_(1J) at frequencies f_(B), f₁+f_(A)−f_(B), f₁+f_(A), and f₁+f_(A)+f_(B) (plus higher-order components, which may be ignored) at point 625. Frequency f_(B) is substantially equal to three intervals (though not identical with frequency f_(A)). Therefore, since frequency f₁+f_(A) is substantially two intervals above f₀, frequency f₁+f_(A)−f_(B) is substantially two intervals below f₀ and frequency f₁+f_(A)+f_(B) is substantially five intervals above f₀.

From point 625, signal Link_(1J) propagates into an amplifier 686. From amplifier 686, signal Link_(1J) propagates through point 624 and into a low-center bandpass filter 688. Like bandpass filter 666, bandpass filter 688 has a passband slightly greater than two intervals and centered on a frequency two intervals below frequency f₀. Since frequency f₁+f_(A)−f_(B) is one interval below frequency f_(O), while frequency f_(B) is more than four intervals below f₀ and frequencies f₁+f_(A) and f₁+f_(A)+f_(B) are above frequency f₀, frequency f₁+f_(A)−f_(B) is passed and frequencies f_(B), f₁+f_(A), and f₁+f_(A)+f_(B) are attenuated.

From bandpass filter 688, signal Link_(1J) propagates through point 625 and into an amplifier 690. Amplifiers 680, 686, and 690 also serve as unidirectional isolators. From amplifier 690, signal Link_(1J) propagates through point 628 and into combiner/splitter 678. The circuit from combiner/splitter 678, through amplifier 680, through filter 682, through mixer 684, through amplifier 686, through filter 688, through amplifier 690, and back to combiner/splitter 678 also forms a loop inhibited from oscillation by filters 682 and 688 in conjunction with mixer 684.

From combiner-splitter 678, signal Link_(1J) propagates through point 629 (coincident with point 620), through combiner/splitter 676, through point 630, and into bandpass filter 692. Like bandpass filter 658, bandpass filter 692 has a passband slightly greater than two intervals and centered on frequency f₀. Since frequency f₁+f_(A)−f_(B) is one interval below frequency f_(O), frequency f₁+f_(A)−f_(B) is passed.

From bandpass filter 692, signal Link_(1J) propagates through point 631, over Link2 signal path 240 through high-frequency propagation medium 202 to second probe 200, through point 632, and into bidirectional hybrid mixer 642. Mixer 642 mixes signal Link_(1J) at frequency f₁+f_(A)−f_(B) from point 632 with signal _(S)IF₃₃ at frequency f_(IF) from point 633 to produce signal Link_(1J) at frequencies f₁+f_(A)−f_(B), f_(IF)−f₁+f_(A)−f_(B), f_(IF), and f_(IF)+f₁+f_(A)−f_(B) (plus higher-order components, which may be ignored) at point 634. Frequencies f₁+f_(A)−f_(B), f_(IF)−f₁+f_(A)−f_(B), and f_(IF) are filtered by receiver 400 (FIG. 5).

From point 634, signal Link_(1J) propagates through combiner/splitter 238 and to point 631. Since frequency f_(IF)+f₁+f_(A)−f_(B) is substantially equal to frequency f_(L1), signal Link_(1J) has successfully been amplified in its passage from first probe 100 to second probe 200.

Signal Link2j propagates from second probe 200, through reciprocal-amplification circuit 600, and into first probe 100 in an analogous manner. The following abbreviated discussion traces signal Link_(2J) through reciprocal-amplification circuit 600 in conjunction with Table 7. Table 7 depicts the frequencies of signal Link_(2J) at each of points 602–634 and 636 throughout reciprocal amplification circuit 600 (points 601 and 635 are not used for signal Link_(2J)).

TABLE 7 Points Frequencies 636 634 f_(L2) 633 f_(IF) 632 631 f_(IF) − f_(L2), f_(IF), f_(L2), f_(IF) + f_(L2) 630 620 621 622 623 f₂, i.e., f_(IF) − f_(L) ₂ 624 f_(B) 625 626 f_(B), f₂ − f_(B), f₂, f₂ + f_(B) 627 628 629 619 618 f₂ − f_(B) 608 609 610 611 612 f_(A) 613 614 f_(A), f₂ − f_(B) − f_(A), f₂ − f_(B), f₂ − f_(B) + f_(A) 615 616 617 607 606 f₂ f_(B) + f_(A) 605 604 f_(IF) 603 602 f₂ − f_(B) + f_(A), f_(IF) − f₂ − f_(B) + f_(A), f_(IF), f_(IF) + f₂ − f_(B) + f_(A) Result f_(IF) + f₂ − f_(B) + f_(A)≈f_(L2)

Signal Link_(2J) originates in second probe 200 with frequency f_(L2). Signal Link_(2J) propagates from point 636, through combiner/splitter 238, through point 634, and into bidirectional hybrid mixer 642. Mixer 642 mixes signal Link_(2J) with signal _(S)IF₃₃ at frequency f_(IF) from point 633 to produce signal Link_(2J) at frequencies f_(L2), f_(IF)−f_(L2), f_(IF), and f_(IF)+f_(L2) at point 632. From point 632, signal Link_(2J) propagates over Link2 signal path 240 through high-frequency propagation medium 202 and into high-low subcircuit 652 at point 631.

From point 631, signal Link_(2J) propagates through bandpass filter 692 to point 630. Bandpass filter 692 has a passband slightly greater than two intervals and centered on frequency f₀. Since f₂=f_(1H)−f_(L2), bandpass filter 692 passes signal Link_(2J) at frequency f₂ while attenuating the other components.

From point 630, signal Link_(2J) propagates through combiner/splitter 676, through point 620, through combiner/splitter 678, through point 621, and into amplifier 680. From amplifier 680, signal Link_(2J) propagates through point 622 and into high-center bandpass filter 682. Bandpass filter 682 has a passband slightly greater than two intervals and centered on a frequency two intervals above frequency f₀. Since frequency f₂ is one interval above frequency f₀, bandpass filter 682 passes signal Link_(2J).

From bandpass filter 682, signal Link_(2J) propagates through point 623 and into mixer 684. Mixer 684 mixes signal Link_(2J) at frequency f₂ with a signal (not shown) at frequency f_(B) from point 624 to produce signal Link_(2J) at frequencies f_(B), f₂−f_(B), f₂, and f₂+f_(B) at point 625. Since frequency f₂ is substantially one interval above f₀, frequency f₂−f_(B) is substantially two intervals below f₀ and frequency f₂+f_(B) is substantially four intervals above f₀.

From point 625, signal Link_(2J) propagates into amplifier 686. From amplifier 686, signal Link_(2J) propagates through point 626 and into low-center bandpass filter 688. Since frequency f₂−f_(B) is two intervals below frequency f₀, frequency f₂−f_(B) is passed while other components are attenuated.

From filter 688, signal Link_(2J) propagates through point 627 and into amplifier 690. From amplifier 690, signal Link_(2J) propagates through point 628, through combiner/splitter 678, through point 629, through combiner splitter 676, through point 619, through propagation medium 654 via LinkP signal path 656, into low-high subcircuit 650, through point 618, through combiner/splitter 660, through point 608, through combiner/splitter 662, through point 609, and into amplifier 664. From amplifier 664, signal Link_(2J) propagates through point 610 and into low-center bandpass filter 666. Since frequency f₂−f_(B) is two intervals below frequency f₀, frequency f₂−f_(B) is passed.

From bandpass filter 666, signal Link_(2J) propagates through point 611 and into mixer 668. Within mixer 668, signal Link_(2J) at frequency f₂−f_(B) is mixed with a signal (not shown) at frequency f_(A) from point 612 to produce signal Link_(2J) at frequencies f_(A), f₂−f_(B)−f_(A), f₂−f_(B), and f₂−f_(B)+f_(A) at point 613. Since frequency f₂−f_(B) is substantially two intervals below f₀, frequency f₂−f_(B)−f_(A) is substantially five intervals below f₀ and frequency f₂−f_(B)+f_(A) is substantially one interval above f₀.

From point 613, signal Link_(2J) propagates through amplifier 670, through point 614, and into high-center bandpass filter 672. Since frequency f₂−f_(B)+f_(A) is one interval above frequency f_(O), frequency f₂−f_(B)+f_(A) is passed and other components are attenuated.

From bandpass filter 672, signal Link_(2J) propagates through point 615, through amplifier 674, through point 616, through combiner/splitter 662, through point 617, through combiner/splitter 660, through point 607, and into bandpass filter 658. Since frequency f₂−f_(B)+f_(A) is one interval above frequency f_(O), frequency f₂−f_(B)+f_(A) is passed.

From bandpass filter 658, signal Link_(2J) propagates through point 606, over Link1 signal path 140 through high-frequency propagation medium 102 to first probe 100, through point 605, and into bidirectional hybrid mixer 640. Mixer 640 mixes signal Link_(2J) at frequency f₂−f_(B)+f_(A) with signal _(S)IF₃₃ at frequency f_(IF) from point 604 to produce signal Link_(2J) at frequencies f₂−f_(B)+f_(A), f_(IF)−f₂−f_(B)+f_(A), f_(IF), and f_(IF)+f₂−f_(B)+f_(A) at point 603. Frequencies f₂−f_(B)+f_(A), f_(IF)−f₂−f_(B)+f_(A), and f_(IF) are filtered by receiver 400 (FIG. 5).

From point 603, signal Link_(2J) propagates through combiner/splitter 138 and to point 602. Since frequency f_(IF)+f₂−f_(B)+f_(A) is substantially equal to frequency f_(L2), signal Link_(2J) has successfully been amplified in its passage from second probe 200 to first probe 100.

It should be noted that reciprocal amplification circuit 600 will desirably provide sufficient amplification to overcome losses in the Link signal path between probes 100 and 200. those skilled in the art will appreciate that propagation medium 654 may be any required propagation medium, e.g., an RF transmission link, earth-conductance link, etc. The use of other media and/or other schema than those presented herein does not depart from the spirit of the present invention.

FIG. 13 shows a simplified schematic diagram of a circular mixer circuit 700 for TES network analyzer 22 in accordance with an alternative preferred embodiment of the present invention. The following discussion refers to FIG. 13.

The reciprocal link leg (i.e., Link signal path 140/240 through Link channel 401 (FIGS. 3, 4, and 5)) provides the necessary amplitude and phase references between probes 100 and 200 (FIGS. 3 and 4) in TES network analyzer 22 of FIGS. 1–5 and 7. However, other embodiments also allow amplitude and phase references to be established between probes 100 and 200.

One such embodiment is the use of circular mixer circuit 700. In this approach, an accurate phase reference between probes 100 and 200 is established by measuring the phase at each port relative to the phase in analysis unit 28. This is accomplished by measuring the time-varying transmission coefficients for RefF1 and RefF2 signal paths 114 and 214, i.e., by measuring T^(RefF1) and T^(RefF2) instead of T^(Link).

As was the case when measuring transmission coefficient T^(Link) for Link signal path 140/240, transmission coefficients T^(RefF1) and T^(RefF2) associated with RefF1 and RefF2 signal paths 114 and 214, respectively, may be measured by simultaneously passing two signals, of different frequencies, in opposite directions through each of RefF1 and RefF2 signal paths 114 and 214. To do this, RefF1 and RefF2 signal paths 114 and 214 will contain reciprocal components. Therefore, RefF1 and RefF2 signal paths 114 and 214 may not include unidirectional components, i.e., amplifiers and isolators. If necessary, duplicate reciprocal-amplification circuits 600, discussed hereinbefore, may be used in RefF1 and RefF2 signal paths 114 and 214.

In circular mixer circuit 700 of FIG. 13, fixed-reference signal RefF₁₃ is generated in analysis unit 28 (i.e., in transmitter 300), passes through directional coupler 702, and propagates through high-frequency propagation medium 102 over RefF1 signal path 114 into first probe 100, where it becomes signal RefF₁₁. Signal RefF₁₁ is split by combiner/splitter 704 into fixed-common signal ComfF₁₁ and a circular-common signal ComC₁₁. Mixing circuit 706 then mixes signal ComC₁₁ with a local-oscillator offset signal _(L)ComC₁₁ (generated by local-oscillator module 354 (not shown)) to produce an offset circular-common signal _(O)ComC₁₁ (see Equations (3) and (4)). Signal _(O)ComC₁₁ passes through combiner/splitter 704 to become a circular signal Circ₁₁. Signal Circ₁₁ propagates through high-frequency propagation medium 102 over RefF1 signal path 114 into analysis unit 28, where it becomes a signal Circ₁₃. Directional coupler 702 extracts signal Circ₁₃ from RefF1 signal path 114 as an extracted circular signal _(X)Circ₁₃. Signal _(X)Circ₁₃ passes through combiner 708, where it is combined with a signal _(X)Circ₂₃ from probe 200 to become a composite circular signal _(X)Circ₃₃ in receiver 400.

Similarly, fixed-reference signal RefF₂₃ is generated in analysis unit 28 (i.e., in transmitter 300), passes through directional coupler 710, and propagates through high-frequency propagation medium 202 over RefF2 signal path 214 into second probe 200, where it becomes signal RefF₂₂. Signal RefF₂₂ is split by combiner/splitter 712 into fixed-common signal ComF₂₂ and a circular-common signal ComC₂₂. Mixing circuit 714 then mixes signal ComC₂₂ with a local-oscillator offset signal _(L)ComC₂₂ (generated by local-oscillator module 354 (not shown)) to produce an offset circular-common signal _(O)ComC₂₂ (see Equations (3) and (4)). Signal _(O)ComC₂₂ passes through combiner/splitter 712 to become a circular signal Circ₂₂. Signal Circ₂₂ propagates through high-frequency propagation medium 202 over RefF2 signal path 214 into analysis unit 28, where it becomes a signal Circ₂₃. Directional coupler 710 extracts signal Circ₂₃ from RefF2 signal path 214 as an extracted circular signal _(X)Circ₂₃. Signal _(X)Circ₂₃ passes through combiner 708, where it is combined with signal _(X)Circ₁₃ from probe 100 to become composite circular signal _(X)Circ₃₃ in receiver 400.

Since Link signal path 140/240 no longer exists, four-way incident-signal couplers 142 and 242 (FIGS. 3 and 4) have been replaced with three-way couplers 716 and 718, respectively.

It may be seen from the above discussion that RefF1 and RefF2 signal paths 114 and 214 are bidirectional signal paths, with signals RefF_(1J) and RefF_(2J) propagating from analysis unit 28 to probes 100 and 200, while signals Circ_(1J) and Circ_(2J) substantially simultaneously propagate from probes 100 and 200 to analysis unit 28, respectively.

Each signal therefore passes through the S-matrix [S^(RefF1)] or [S^(RefF2)]. Following the procedures that were used to create Tables 3 and 4, we find that signals Circ_(1J) and Circ_(2J) may be written as _(M) ComC ₁₃ =A _(S)(T ^(RefF1))² P _(MComC13) exp(−j∠A _(LComC11))G ^(Circ), and  (41) _(M) ComC ₂₃ =A _(S)(T ^(RefF2))² P _(MComC23) exp(−j∠A _(LComC22))G ^(Circ).  (42)

Taking the square root of the ratio of these two signals gives,

$\begin{matrix} {\frac{T^{RefF2}}{T^{RefF1}} = {\sqrt{\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{P_{MComC13}}{P_{MComC23}}\frac{\exp\left( {\angle\; A_{LComC22}} \right)}{\exp\left( {\angle\; A_{LComC11}} \right)}}.}} & (43) \end{matrix}$

Signal normalization techniques can now be used to compute the forward and reverse, measured transmission coefficients of DUT 26 using

$\begin{matrix} \begin{matrix} {S_{21\; M} = {\frac{T^{RefF2}}{T^{RefF1}}\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{\exp\left\lbrack {j\left( {{\angle\; A_{LComF11}} + {\angle\; A_{LDat22}} + {\angle\; A_{LComV22}}} \right)} \right\rbrack}{\exp\left\lbrack {j\left( {{\angle\; A_{LInc11}} + {\angle\; A_{LComD22}} + {\angle\; A_{LComF22}}} \right)} \right\rbrack}}} \\ {{= {{\overset{\sim}{S}}_{21\; M}\frac{P_{MComF13}}{P_{MInc13}}\frac{P_{MDat23}}{P_{MComD23}}\frac{P_{MComV23}}{P_{MComF23}}}},{and}} \end{matrix} & (44) \\ \begin{matrix} {\;{S_{12M} = {\frac{T^{RefF1}}{T^{RefF2}}\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{{}_{}^{}{}_{}^{}}{{}_{}^{}{}_{}^{}}\frac{\exp\left\lbrack {j\left( {{\angle\; A_{LComF22}} + {\angle\; A_{LDat11}} + {\angle\; A_{LComV11}}} \right)} \right\rbrack}{\exp\left\lbrack {j\left( {{\angle\; A_{LInc22}} + {\angle\; A_{LComD11}} + {\angle\; A_{LComF11}}} \right)} \right\rbrack}}}} \\ {= {{\overset{\sim}{S}}_{12\; M}\frac{P_{MComF23}}{P_{MInc23}}\frac{P_{MDat13}}{P_{MComD13}}{\frac{P_{MComV13}}{P_{MComF13}}.}}} \end{matrix} & (45) \end{matrix}$

The reflection coefficients of DUT 26 can still be computed using Equations (13) and (16).

It is interesting to compare Equations (43), (44), and (45) with Equations (14), (15), and (18). These represent various options for establishing the amplitude and phase references between probes 100 and 200. All such options, however, depend on having a reciprocal signal flow path. It has been demonstrated herein that Link signal path 140/240 and RefF1 and RefF2 signal paths 114 and 214 may be used for this purpose. Those skilled in the art will appreciate that it is also possible to employ either RefV1 and RefV2 signal paths 118 and 218 or Tx1 and Tx2 signal paths 112 and 212 for this purpose. The use of other signal paths to effect amplitude and phase references between probes 100 and 200 does not depart from the spirit of the present invention.

FIG. 14 shows a simplified schematic diagram of a mix/demix circuit in accordance with an alternative embodiment of the present invention. The following discussion refers to FIG. 14.

As shown in mix/demix circuit 800 of FIG. 14, the TES technique can be adapted for use with a standard (i.e., prior-art) network analyzer 802. Note that circuit 800 only allows for the measurement of the transmission coefficient in a single direction through DUT 26. Circuit 800 also assumes that accurate amplitude and phase references exist on ports 30 and 32 of DUT 26, i.e., no link leg is required. Only time-varying transmitter errors [S^(Tx)] and receiver errors [S^(Rx)] are simulated. Those skilled in the art will appreciate that circuit 800 is exemplary only, and that it would be straightforward to expand circuit 800 to a full two-port network-analyzer measurement system. Such expansion of circuit 800 does not depart from the spirit of the present invention.

Signal source 24 produces source signal _(S)Tx₃₃ at predetermined frequency f_(S). Signal _(S)Tx₃₃ is transmitted through time-varying S-matrix [S^(Tx)]. S-matrix [S^(Tx)] models the time-varying errors associated with the transmitter circuitry. A splitter 804 splits signal _(S)Tx₃₃ into two paths. One of these paths flows through DUT 26 and the other flows through two frequency-offset mixers 806 and 808. The output signal from mixer 808 appears at the four sum and difference frequencies f _(S) ±f _(L1) ±f _(L2),  (46) where f_(L1) and f_(L2) denote the frequencies for two frequency-offset local-oscillator sources 810 and 812. Since these signals are at different frequencies than signal _(S)Tx₃₃ that flows through DUT 26, combiner 814 can combine the signals in the two paths. The combined signals then flow through a pad 816 and an amplifier 818 before reaching time-varying S-matrix [S^(Rx)] that models the time-varying receiver errors. Pad 816 and amplifier 818 are included to isolate mixers 806 and 808 and DUT 26 from time-varying S-matrix [S^(Rx)].

After passing through S-matrix [S^(Rx)] (i.e., a receiver cable), the signals must be detected. In TES network analyzer 22, the time-domain data and normalization signals are both detected simultaneously by analog-to-digital converters 493–496 (FIG. 5). FFT routine 510 in processor 500 (FIG. 7) is then employed to transform the time-domain signals into the frequency domain, where the individual data and normalization signals are extracted from the appropriate frequency bins. However, a receiver (not shown) of standard network analyzer 802 is configured to only detect a single CW wave. If a narrow enough resolution bandwidth (i.e., RBW<<min(f_(L1),f_(L2))) is selected for the IF filter in standard network analyzer 802, then the signal that is detected in an A-channel 820 of network analyzer 802 will provide information about the transmission characteristics of DUT 26. Likewise, after passing through two more mixers 822 and 824, an R-channel 826 will provide information about the incident signal. The desired transmission coefficient can then be obtained via signal normalization, i.e.,

$\begin{matrix} {S_{21} = {\frac{{Dat}_{A}}{{Dat}_{R}}.}} & (47) \end{matrix}$

There are a number of issues that should be noted. First, the signal normalization removes time-varying errors associated with source-drift, transmitter circuitry drift, and the receiver circuitry drift up to a splitter 828. However, any drift errors associated with the receiver circuitry of network analyzer 802 are not removed, since the data and normalization signals are detected in separate receiver channels. It should also be noted that mixers 806, 808, 822, and 824 are used to reduce RF feedthrough to a level that doesn't cause detectable errors in the measured signals.

The signals from local-oscillator sources 810 and 812 are directed to mixer-pairs 806/824 and 808/822, respectively, by splitters 830 and 832, and 808 and 822, respectively. Pads 816, 834, 836, 838, 840, and 842 are used to provide isolation and impedance matching as required.

This technique, which is referred to as a mix/demix technique, may be extended to full two-port DUT measurements. However, it would require numerous receiver channels in order to simultaneously measure all S-Parameters. There are many applications where mix/demix circuit 800 may be cost-effectively used with standard network analyzer 802. For example, when the goal is to measure a single quantity, such as a direct measurement from one transducer to another, there are sufficient receiver channels in a standard network analyzer for the complete measurement. If the transducers are remote from the network analyzer, virtually all of the error due to time-varying effects may be in the link between the transducers and network analyzer 802. In this case the time-varying errors in the separate data-acquisition channels may be negligible. The mix/demix implementation of the TES technique normalization can then be a very powerful and cost-effective approach.

FIG. 15 shows a simplified schematic diagram of an intermodulation distortion checking and removal circuit in accordance with a preferred embodiment of the present invention. The following discussion refers to FIG. 15.

As discussed hereinbefore, care must be taken to avoid intermodulation distortion (IMD) since multiple signals existing at different frequencies simultaneously pass through various nonlinear elements, e.g., mixers and amplifiers. An IMD program may be used to minimize the effects of IMD by properly choosing the frequencies of the local-oscillator signals (see Table 1) for the frequency-offset mixers during the design of TES network analyzer 22. Unfortunately, IMD programs only account for the effects of the ideal fundamental frequencies and the associated harmonics. They do not account for the effects of spurious signals that are either generated in signal measurement system 20 or caused by electromagnetic interference from strong external signals. Furthermore, it is only possible to design system 20 to avoid IMD for intermodulation products that are less than some finite order and it is difficult to predict the actual IMD level. Finally, if DUT 26 is nonlinear, IMD during signal measurement may lead to errors.

Since IMD effects are unpredictable and may be influenced by the measurement setup and external environmental factors (e.g., power-line noise and electromagnetic interference), a test that indicates the level of IMD in a given measurement is highly desirable. This may be accomplished by changing the level of the signal input into the nonlinear components. The measured results for the two levels may then be compared to indicate the level of IMD and nonlinearity. For example, if the signal levels input into an amplifier are changed by 3 dB, then the desired output signal levels, which should be linearly related to the input signal levels, should also change by 3 dB. A deviation from ideal linear behavior indicates that there is either IMD (e.g., second-order and third-order intermodulation products should change by approximately 6 dB and 9 dB, respectively) or signal compression. IMD test and reduction circuit 900 (FIG. 15) may be used to test for IMD and nonlinearities. In order to simplify FIG. 15, only those components of first probe 100 (see FIG. 3) specifically a part of circuit 900 are shown). Those skilled in the art will appreciate that while first probe 100 is exemplified, a substantially identical circuit may be incorporated into second probe 200. The use of circuit 900 in either or both probes 100 and/or 200 does not depart from the spirit of the present invention.

When switches 902, 904, 906, and 908, are set to “through” positions (i.e., the positions shown), raw data signal _(R)Dat₁₁ is divided in a hybrid splitter 910 into signals _(X)Dat₁₁ and _(X)IMod₁₁. Signals _(X)Dat₁₁ and _(X)IMod₁₁ flow into mixing circuits 128 and 912, where they are mixed with frequency-offset local-oscillator signals _(L)Dat₁₁ and _(L)IMod₁₁, respectively to produce offset signals _(O)Dat₁₁ and _(O)IMod₁₁. Local-oscillator signals _(L)Dat₁₁ and _(L)IMod₁₁ exist at different frequencies, so offset signals _(O)Dat₁₁ and _(O)IMod₁₁ are offset from transmission signal Tx₁₁ (i.e., source signal _(S)Tx₃₃ (FIG. 2) by different amounts. Offsetting signals _(O)Dat₁₁ and _(O)IMod₁₁ by different amount allows them to be recombined by combiner 914 into composite data signal Dat₁₁.

If pads 129 and 913 (input pads 36 (FIG. 6) of mixing circuits 128 and 912) have different values, e.g., were mixing circuit 128 to attenuate signal _(X)Dat₁₁ 3 dB less than mixing circuit 912 attenuates signal _(X)IMod₁₁, then the levels of signals _(X)Dat₁₁ and _(X)IMod₁₁, the signal levels that pass through the nonlinear components, i.e., mixers 127 and 911 (mixers 42 (FIG. 6) of mixing circuits 128 and 912), optional amplifier 156, and the components of Dat1 channel 404 of receiver 400 (FIG. 5) also differ by 3 dB. These two signals, being simultaneously measured, may be used to efficiently test for IMD and nonlinearities. If pad 913 were switched in, instead of being fixed in circuit 900, then potential switch repeatability errors would necessitate system calibration after each measurement. Also, the use of a switch would necessitate doubling the number of measurements, leading to lower productivity and a significant increase in measurement cost.

It is also possible to check for IMD and nonlinearities by changing the power level of transmission signal Tx₁₁. However, this would change the power level seen by DUT 26 and would require twice as much measurement time.

In the alternative exemplary embodiment of FIG. 15, IMD testing is shown applied to the data legs, as this is the most-likely location to suffer IMD effects. Those skilled in the art will appreciate that the schema of circuit 900 may also be used in the incident legs if desirable. The use of this schema in other than the data legs does not depart from the spirit of the present invention.

On the rare occasion when such IMD and nonlinearities are detected, a method is presented for removing those errors. This is accomplished by placing switches 902, 904, 906, and 908 in “attenuated” positions (opposite the positions shown). This places attenuators 916 and 918 into the paths of data and common signals _(R)Dat₁₁ and ComD₁₁. Attenuators 916 and 918 are placed into both signal paths to allow the measured data, intermodulation, and data-common signals _(M)Dat₁₃, _(M)IMod₁₃, and _(M)ComD₁₃ within composite data signal Dat₁₁ to be maintained within ±15 dB of one another as described hereinbefore.

Placing switches 902, 904, 906, and 908 in the “attenuated” positions lowers the levels of signals _(R)Dat₁₁ and ComD₁₁, and reduces the effects of IMD and nonlinearities. Circuit 900 shows single pads 916 and 918 of up to 20 dB placed into the path of each of signals _(R)Dat₁₁ and ComD₁₁, respectively. Those skilled in the art, however, will appreciate that switches 902, 904, 906, and 908 and attenuators 916 and 918 may be replaced by step switches with multiple attenuators, or by variable attenuators, if more flexibility in attenuation is desired, without departing from the spirit of the present invention.

The effects associated with switching between the “through” and “attenuated” positions are not removed by normalization, since multiple signals do not flow through these components. Thus, a new static calibration may be required after switching.

In summary, the present invention teaches a real-time error-suppression method and apparatus therefor. The combination of static, dynamic, phase-noise, and random error suppression incorporated into a Thorough Error Suppression (TES) technique leads to a significant improvement in signal measurement accuracy. The TES technique also makes it possible to use much-lower-cost components in the construction of measuring systems. Examples of these savings are lower cost and simpler components in the receiver, the use of conventional over high-cost specialized cables connecting all devices and ports, and conventional off-the-shelf signal-conditioning circuitry. This leads to a significant economic benefit by calibrating less expensive components, rather than building very expensive, precisely matched components.

Although the preferred embodiments of the invention have been illustrated and described in detail, it will be readily apparent to those skilled in the art that various modifications may be made therein without departing from the spirit of the invention or from the scope of the appended claims. 

1. A method of noise suppression in a signal measurement system, said method comprising: a) providing a first signal; b) deriving a second signal from said first signal, wherein said second signal is an analog signal; c) deriving a third signal from said first signal, wherein said third signal is an analog signal; d) analyzing said third signal to determine the presence of noise within said first signal; and e) compensating, when said analyzing activity d) determines said noise is present within said first signal, said second signal in response to said noise within said first signal.
 2. A method as claimed in claim 1 wherein said deriving activities b) and c) are executed substantially concurrently.
 3. A method as claimed in claim 1 wherein: said method additionally comprises, prior to said compensating activity e), converting said second signal into a fourth signal, wherein said fourth signal is a digital signal; said method additionally comprises, prior to said analyzing activity d), converting said third signal into a fifth signal, wherein said fifth signal is a digital signal; said analyzing activity d) analyzes said fifth signal; and said compensating activity e) compensates said fourth signal.
 4. A method as claimed in claim 3 wherein said analyzing activity d) comprises: transforming said fifth signal from the time domain to the frequency domain to produce a sixth signal distributed into a first plurality of frequency bins; and categorizing each of said frequency bins as a noise bin containing a noise energy peculiar thereto.
 5. A method as claimed in claim 4 wherein said compensating activity e) comprises: transforming said fourth signal from the time domain to the frequency domain to produce a seventh signal distributed into a second plurality of frequency bins; correlating each frequency bin in said second plurality with a corresponding frequency bin in said first plurality; and subtracting said noise energy of each frequency bin in said first plurality from said correlated frequency bin in said second plurality.
 6. A method as claimed in claim 1 wherein said deriving activity c) comprises: extracting a fourth signal from said first signal; providing a fifth signal; mixing said fourth signal and said fifth signal to produce a sixth signal; and down-converting said sixth signal to produce said third signal.
 7. A method as claimed in claim 1 wherein: said method additionally comprises calibrating said third signal to determine invariant phase perturbations therefor; and said compensating activity e) comprises subtracting said invariant phase perturbations from said second signal.
 8. A method as claimed in claim 1 wherein said deriving activity b) comprises: extracting a fourth signal from said first signal; providing a fifth signal; mixing said fourth signal and said fifth signal to produce a sixth signal; and down-converting said sixth signal to produce said second signal.
 9. A method as claimed in claim 1 wherein: said providing activity a) provides said first signal at a first frequency said deriving activity b) derives said second signal at a second frequency; said deriving activity c) derives said third signal at a third frequency; and said method additionally comprises establishing a frequency band, all frequencies of which are within ten percent of a center frequency thereof, wherein each of said first, second, and third frequencies is a different frequency within said frequency band.
 10. A method as claimed in claim 1 wherein: said method additionally comprises producing a fourth signal; said method additionally comprises injecting said fourth signal into said first signal to produce a fifth signal, said fifth signal being a composite of said first signal and said fourth signal; said deriving activity c) derives said third signal from said fifth signal; and said method additionally comprises determining invariant phase perturbations of said first signal from said third signal.
 11. A method as claimed in claim 1 additionally comprising: coupling a first probe to a first port of a device-under-test; coupling a second probe to a second port of said device-under-test; propagating said first signal from said first probe into said first port of said device-under-test to said second port of said device-undertest and into said second probe; deriving, within said first probe, a fourth signal from said first signal; deriving, within said first probe, a fifth signal from said first signal; deriving, within said second probe, a sixth signal from said first signal; deriving, within said second probe, a seventh signal from said first signal; reciprocally propagating said fifth and seventh signals from said first and second probes to said second and first probes, respectively, over a common signal path; combining, within said first probe, said fourth signal and said seventh signal to produce said second signal; and combining, within said second probe, said sixth signal and said fifth signal to produce an eighth signal.
 12. A method as claimed in claim 11 additionally comprising: interposing a reciprocal amplification circuit into said common signal path between said first and second probes; amplifying said fifth signal during said propagation from said first probe to said second probe; and amplifying said seventh signal during said propagation from said second probe to said first probe.
 13. A method as claimed in claim 1 wherein said noise is phase noise.
 14. A method of calibrating a noise channel in a signal measurement system, said method comprising: providing a first signal; producing a second signal; injecting said second signal into said first signal to produce a third signal, said third signal being a composite of said first signal and said second signal; deriving a fourth signal from said third signal; and determining invariant perturbations of said first signal from said fourth signal.
 15. A method as claimed in claim 14 wherein said producing activity comprises: extracting a fifth signal from said first signal; providing a sixth signal; and mixing said fifth and sixth signals to produce said second signal.
 16. A method as claimed in claim 15 wherein said deriving activity comprises: extracting a seventh signal from said third signal; providing an eighth signal; mixing said seventh and eighth signals to produce a ninth signal; and down-converting said ninth signal to produce said fourth signal.
 17. A method as claimed in claim 14 wherein, prior to said determining activity, said fourth signal is an analog signal, and wherein said determining activity comprises: converting said fourth signal into a fifth signal, said fifth signal being a digital signal in the time domain; transforming said fifth signal from the time domain to the frequency domain to produce a sixth signal distributed into a plurality of frequency bins; determining which of said frequency bins contain energy of said perturbations, each of said bins being a noise bin containing a noise energy peculiar thereto; and categorizing each of said frequency bins as a noise bin containing a noise energy peculiar thereto.
 18. A method as claimed in claim 14 wherein said perturbations comprise noise.
 19. A method as claimed in claim 18 wherein said noise is phase noise.
 20. A method of noise suppression in a signal measurement system, said method comprising: providing a first signal; deriving a second signal from said first signal, wherein said second signal is an analog signal; deriving a third signal from said first signal substantially concurrently with said second signal, wherein said third signal is an analog signal; converting said second signal into a fourth signal, wherein said fourth signal is a digital signal; converting said third signal into a fifth signal, wherein said fifth signal is a digital signal; analyzing said fifth signal to determine the presence of phase noise within said first signal; and compensating, when said analyzing activity d) determines said phase noise is present within said first signal, said fourth signal in response to said noise within said first signal. 